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HAUPTAUFSÄTZE 


Die Bestimmung des Sicherheitsgrades einfach statisch 
unbestimmter Fachwerke. 


Von Ernst Melan in Wien. 


1. Der Sicherheitsgrad eines statisch unbestimmten Fachwerkes. Läßt man alle äußeren 
Kräfte, welche auf eine Tragkonstruktion einwirken, proportional wachsen, so wird bei einem 
bestimmten Vielfachen der gegebenen Belastung das System unbrauchbar, das heißt, es werden 
Formänderungen desselben ohne Aufwand äußerer Arbeit möglich. Man nennt den Faktor, 
mit welchem die tatsächlich einwirkende Belastung zu multiplizieren ist, damit dies eintritt, 
den Sicherheitsgrad der betreffenden Konstruktion. 

Bei einem statisch bestimmten System ist die Bestimmung des Sicherheitsgrades sehr 
einfach. Speziell bei einem statisch bestimmten Fachwerk ist der Sicherheitsgrad der ganzen 
Konstruktion gleich dem des schwächsten Stabes. Denn ein statisch bestimmtes Fachwerk 
hat gerade nur so viele Stäbe als unbedingt notwendig sind, um die gegenseitige Lage der 
Knoten festzulegen. Ist also ein einziger Stab unbrauchbar geworden, d. h. sind Änderungen 
seiner Länge und damit der Entfernung der beiden durch ihn verbundenen Knoten möglich, 
ohne daß die Stabkraft wächst, so ist das ganze Fachwerk unbrauchbar geworden. 

Wesentlich anders liegen die Verhältnisse bei einem statisch unbestimmten Fachwerk. 


Damit ein r-fach statisch unbestimmtes Fachwerk unbrauchbar wird, ist es im allgemeinen 


notwendig, daß r +1 Stäbe unbrauchbar werden. Denn ein solches System hat ja r Stäbe 
mehr, als zur Festlegung der gegenseitigen Lage der Knoten notwendig sind. Es müssen 
daher, keine Häufungsstellen der Stäbe vorausgesetzt, r-+1 Stäbe unbrauchbar werden, da- 
mit der Rest des Fachwerkes unbrauchbar wird. 

Man pflegt in der Praxis dem Erfordernis eines bestimmten Sicherheitsgrades in der 
Weise Rechnung zu tragen, daß man eine „zulässige* Spannung vorschreibt, welche 


einen Bruchteil, etwa „, jener beträgt, bei welcher der betreffende Stab unbrauchbar wird. 


Bei einem statisch bestimmten Fachwerk ist hiermit auch der Sicherheitsgrad des ganzen 
Systemes, nämlich gleich n, gegeben, wenn die zulässige Spannung auch nur in einem ein- 
zigen Stabe erreicht wird. Bei einem r-fach statisch unbestimmten Fachwerk ist es aber 


notwendig, daß die zulässige Spannung im allgemeinen in mindestens r +1 Stäben erreicht 


= 7; 
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wird. Denn dann erst werden bei Anwachsen der Belastung auf den »-fachen Wert r+1 
Stäbe und damit das System unbrauchbar, der Sicherheitsgrad somit gleich » sein '). 

Ist die Bedingung, daß in r +1 Stäben die zulässige Spannung erreicht wird, wie dies 
häufig der Fall ist, nieht erfüllt, so kann aus der zulässigen Spannung nicht ohne weiteres 
auf den Sieherheitsgrad der Konstruktion geschlossen werden. Es ist vielmehr die Aufgabe 
der vorliegenden Untersuehung, den Sicherheitsgrad eines solehen Fachwerkes zu bestimmen, 
wobei wir uns überdies auf einfach statisch unbestimmte Fachwerke beschränken wollen. 

Zunächst sei die Voraussetzung gemacht, daß zwischen der Änderung des Abstandes 
der beiden Knoten eines Stabes und seiner Stabkraft der in Abb. I dargestellte Zusammen- 
hang bestehe. Naeh anfänglieher Proportionalität zwischen Längenänderung und Stabkraft, 
also nach rein elastischem Verhalten des Materiales, folgt nach Erreichen der Traggrenze bei As’ 
ein Bereich, in welehem die Stabkraft den konstanten Wert T besitzt, um schließlich nach Er- 
reichen des Wertes Is” der Längenänderung auf Null zu sinken. Das Verhältnis der Längen- 
änderung an dieser Stelle zu jener an der Traggrenze sei mit g bezeichnet. Für einen Zug- 
stab entspricht unsere Annalıme einem idealplastischen Baustoff, und Abb. I ist bis auf den 
Maßstab mit der Spannungs-Dehnungs-Kurve identisch. Der Wert von g beträgt für Eisen, 
das als Baustoff fast ausschließlich in Frage kommt, ungefähr 200. Für Druckstäbe nehmen 
wir prinzipiell denselben Zusammenhang zwischen Längenänderung und Stabkraft an. Natür- 
lieh wird im allgemeinen mit Rücksicht auf die Kniekung die Tragfähigkeit für Druck, T’, 
kleiner sein wie für Zug. Aber auch der Wert von g ist hier viel kleiner, ja er kann selbst 
gleich 1 werden, wenn der Stab nach Erreiehen der Knieklast sich sofort gänzlich ausschaltet 
und überhaupt die Fähigkeit verloren hat, eine Kraft zu übernehmen. 


45’ 25” / 
RA1T322 
Abb. 1. Abb, 2; 


Nach den eingangs getroffenen Festsetzungen wird ein Stab bei Erreichen der 'Trag- 
erenze unbrauehbar. 


2. Die Stabkräfte eines einfach statisch unbestimmten Fachwerkes nach Erreichen der 
Traggrenze eines Stabes. Die Stabkraft in dem Stabe i eines einfach statisch unbestimmten 
Fachwerkes mit an Stäben kann bekanntlich dureh die Gleichung 


dargestellt werden. Darin bedeuten die 4% die Stabkräfte im Grundsystem, d. h. jenem 
statisch bestimmten Fachwerk, welches man durch Durehschneiden des überzähligen Stabes « 
erhält, S, die Stabkraft in dem überzähligen Stabe a und 0° die Stabkräfte im Grundsystem, 
wenn gesetzt wird. Es ist also speziell AR = 0 und 1. 


Ordnen wir den » Werten S;, A} und @7 die Komponenten der Vektoren 3, Re Qe in 
einem n-dimensionalen Raume zu, so kann das vorstehende Gleichungssystem (1) durch die 
einzige Gleichung 


ersetzt werden, d. h. die Endpunkte aller Vektoren Z und N“ liegen auf einer Geraden q. 
welche die Riehtung von X“ hat. Nachdem bei der Wahl eines anderen Stabes, etwa des 


!, Um Mißverständnisse auszuschließen, sei darauf hingewiesen, daß das Verhalten eines statisch unbestimmten 
Fachwerkes bei wiederholter Belastung und Entlastung nicht Gegenstand dieser Untersuchung ist. Ebenso ist auf 
Gleichgewichtszustände, welehe nach dem Erreiehen der Traggrenze in r+ 1 Stäben 


auf Grund unserer Annahmen 
allerdings nur bei unstetiger Änderung der Verschiebungen - 


möglich sind, nieht Rüeksieht genommen. 
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Stabes b, offenbar dieselben Stabkräfte und daher auch derselbe Vektor Z sich ergeben 


müssen, gilt 


und die sich jetzt ergebende Gerade, die die Riehtung von X’ hat, muß mit der vorerwähnten 
Geraden g identisch sein. Demnach haben alle X dieselbe Riehtung, es ist also 


77 

und auch 
(3; 


für jedes und wenn speziell gesetzt wird, wird wegen Q, = 1: 


und insbesondere 
a 7 


In dem erwähnten »-dimensionalen Raume liegen daher alle © und NR in einer Ebene 7, 
welche durch die Gerade qg und den Ursprung gegeben ist. Die Vektoren R” fallen in die 
Schnittgeraden dieser Ebene mit den Hyperebenen 0, da ja ist. 

Wächst die gegebene Belastung, für welehe Z und R bestimmt wurden, auf den »-fachen 
Wert, so gilt 


denn die Stabkräfte im Grundsystem erhalten jetzt den »-fachen Wert » RY, während die 0% 
unverändert bleiben. An Stelle der Geraden g tritt jetzt die zu ihr parallele .q (vr), welche, 
da » N“ dieselbe Richtung wie NR" hat, ebenfalls in der Ebene „ liegt. Abb. 2 zeigt die 
Ebene » mit den ihr liegenden Vektoren (vr) und »R und der Geraden g(). 

Erweitern wir nun unseren n-dimensionalen Raum durch Hinzunahme einer neuen 
Richtung, in welcher wir die Werte von » auftragen, zu einem n» + 1 dimensionalen, so stellt die 
Gl. (3) in diesem Raume eine Ebene e vor, welche wir uns durch die Projektion der Geraden a 
in » I aus dem Ursprung entstanden denken können. Wir wollen diese Ebene durch zwei in 
ihr liegende Gerade, die durch den Ursprung gehen, darstellen, und können sofort ein solches 
Geradenpaar angeben. Wir brauchen dazu bloß die Endpunkte von RN" und & in v» 1 mit 
dem Ursprung zu verbinden und erhalten so die Geraden r* und j. Anstatt der Geraden | 
können wir aber auch die Gerade x” verwenden, die durch die Verbindung des Endpunktes 
von ®”’ mit dem Ursprung entsteht. 

Es soll nun der geometrische Ort aller Endpunkte von (vr) in der Ebene & bestimmt 
werden, wenn sich der Wert von » ändert. Zunächst sei bemerkt, daß in praktisch vor- 
kommenden Fällen das Fachwerk stets so dimensioniert ist, daß die Stabkräfte S; bei der 
tatsächlich vorhandenen Belastung immer in jenem Bereiche liegen, wo zwischen Längen- 
änderung und Stabkraft Proportionalität besteht. In diesem Falle berechnet sich die über- 
zählige Stabkraft bekanntlich aus der Gleichung 


„a 
R; F: 
ı 


Su (vr) = (4). 


worin F; und e; Querschnittsfläche und Länge des Stabes ö bedeutet und damit wird 
Somit ergeben die Koordinaten der oben erwähnten Geraden j an der Stelle » die Werte 
von S; (v) =y S; . 
Hat aber etwa der Stab j die Traggrenze erreicht, so gilt die Gl. (4) nicht mehr. An ihre 
Stelle tritt jetzt gemäß der Voraussetzung eines ideal plastischen Materiales die Bedingung 


welche in Verbindung mit Gl. (3) besagt, daß die Endpunkte von © (vr) nunmehr auf dem 
Schnittgebilde der Hyperebene S;(r)=const mit der Ebene e liegen. Dieses Schnittgebilde 
ist die Gerade j;: denn die Hyperebene eines Raumes schneidet sich mit einer Ebene in einer 
Geraden. Hat also ein Stab die Traggrenze erreicht, so tritt an Stelle der Geraden | die 
Gerade j. » kann nun noch so lange weiter wachsen, bis noch ein zweiter Stab die Traggrenze 
erreicht. Ist dies der Stab I, so brauchen wir bloß nachzusehen, wo sich die Gerade j mit I 
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schneidet. Die Gerade | ist die Schnittgerade der Hyperebene 8; (r)= T;=const mit der 
Ebene e und die Abszisse »;; des Schnittpunktes von j und [ ist der Sicherheitsgrad des Systemes. 

Denkt man sich die Schnittgeraden aller Hyperebenen mit der Ebene & konstruiert, 
so wird aus der Ebene e ein Polygon herausgeschnitten, in dessen Inneren der Koordinaten- 
ursprung liegt. Dabei interessiert, wenn man ausschließt, daß die äußeren Kräfte ihren 
Richtungssinn umkehren, nur jener Teil des Polvgones, welcher auf der Seite positiver v liegt. 
Für unseren Zweck ist es ferner ausreichend, wenn wir die Projektion des Polygones in einer 
beliebigen vk- Ebene kennen, denn es kommt ja nur auf die Koordinaten v, ,„ der Polygonecken an. 

Auf Grund dieser Überlegungen können wir mittels eines einfachen zeichnerischen 
Verfahrens zum Ziele kommen. Zu diesem Zwecke stellen wir den n-+1 dimensionalen 
Raum, welcher, wie schon erwähnt, durch die an Werte von S; und die Werte von » gebildet 
wird, mittels seiner » Projektionen in die v i-Ebenen dar. Tragen wir an der Stelle v=1 
jedesmal die Werte von RF und S; oder statt letzterem auch von R” auf, und verbinden die 
Endpunkte mit dem Ursprung, so erhalten wir die Bilder »/ und s;, bzw. r’ der Geraden r“ 
und j oder vr’, wodurch die Ebene & festgelegt ist. Um diese Ebene mit der Hyperebene 
S;(») = eonst zu schneiden, beachten wir, daß die Projektion der Schnittgeraden i in der 
v i-Ebene einfach als Parallele ö; im Abstande T; zur v-Achse erhalten wird. Will man also 
die Projektion des Polygones in der k-Ebene bestimmen, so braucht man bloß die Schnitt- 
punkte von i; mit den beiden Geradenbildern in der » i-Ebene auf die Geradenbilder in der 
vk- Ebene zu projizieren. Die Verbindung dieser beiden Punkte gibt das Bild i, der Geraden i 
in der v#-Ebene.  Verfährt man mit allen Stäben auf diese Weise, so erhält man die 
Projektion des Polygones in der v%-Ebene. Natürlich müssen sich — und das wird in der 
Regel der Fall sein, nicht alle Stäbe mit ihren Schnittgeraden an der Umrandung des 
Polygones beteiligen. 

In Abb. 3 ist außer der vk-Ebene, welche für die Projektion des Polygones gewählt 
wurde, von den anderen Projektionsebenen des »n +1 dimensionalen Raumes nur noch die 
die » i-Ebene darstellt. Die Konstruktion der übrigen Polygonseiten hat ebenso wie für i; 
und zu geschehen. 


Abh. 3 Abb. 4. 


Wir können nun leieht für jeden Wert der S; den Sicherheitsgrad des Fachwerkes 
angeben. Dazu zeiehnen wir uns das Bild s; der Geraden f, indem wir in » 1 den Wert 
von S;. auftragen und mit dem Ursprung verbinden. Die so erhaltene Gerade s; gilt bis zu 
ihrem Schnittpunkt mit einer der Polygonseiten, etwa mit pz. Von hier an ist dann der zu 
dieser Polvgonseite gehörige Stab p an der Traggrenze angelangt. Mit weiter wachsendem v 
geben die Ordinaten von p7; den Wert von S,(») an, bis bei Erreichen der nächsten Polygon- 
ecke, die etwa von p7; mit q7; gebildet wird, nun auch der Stab y an die Traggrenze gelangt 
ist: damit ist das Fachwerk unbrauchbar geworden. Der Sicherheitsgrad des Systemes ist 
also dureh die Abszisse jener Polygonecken gegeben, welche im Sinne wachsender » dem 
Schnittpunkt von s; unmittelbar folgt. 
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3. Die Abhängigkeit des Sicherheitsgrades von den Stabquerschnitten. Zur Lösung einer 
Aufgabe, welche einige praktische Bedeutung besitzt, wollen wir untersuchen, wie sich der 
Vektor © ändert, wenn wir die Fläche F7, eines Stabes ! ändern. Nehmen wir diesen Stab 
vorderhand als den überzähligen Stab an, so gilt, wie schon erwähnt, 


1ei 


wenn wir 


@; 
i 


gesetzt haben. Dabei bedeutet #, r die Spannung an der Traggrenze im Stabe 1. Die 


Summen in a; und b, erstrecken sich über alle Stäbe des Fachwerkes mit Ausnahme des 
Stabes sie sind demnach von unabhängig. db; ist überdies stets positiv, weil 
dasselbe Vorzeichen wie 8; hat. Wenn also zwei zusammengehörige Werte von 5, und 
F,, oder, was dasselbe ist, von S; und T,, also etwa und gegeben sind, 
so kann man, konstantes f, vorausgesetzt, für jedes T; das zugehörige S; aus der Tatsache 
finden, daß sich gemäß der Gl. (6) alle Geraden s; mit den zugehörigen I, d. i. den im 
Abstande T, Parallelen zur v»-Achse auf einer Geraden 0) schneiden. 6) kann natürlich auch 
so bestimmt werden, daß man den Abstand 5, auf der Abszissenachse aufträgt; die Neigung 
da; 
b; 

Nachdem die Geraden 0} als Bild einer Geraden o! in dem n-+ 1 dimensionalen Raum 
aufgefaßt werden kann, kann die Konstruktion, wie auch aus Abb. 4 ersichtlich, in jeder 
beliebigen Projektionsebene durchgeführt werden. 


von 0) ist dann und 0) ist zu s; für FRA >x parallel. 


Es soll noch untersucht werden, wie sich die Lage von der Geraden 0’ mit der Fläche F; 


eines beliebigen Stabes i ändert. Setzt man den Wert von 2 R Q' Fr. aus der Gleichung 
J 


2,0; 


Se — in Gl. (6a) ein, so ergibt sich für b, der Wert 
F; 
e,tı 
© 
Ss, 20; 
Analog wird 
h;: 
0; 
J 
S; F': 
J 


wobei sich die Summen in den Nennern beider Ausdrücke über alle Stäbe erstrecken. 
Nachdem nach Gl. (2a) 


0 
ist, folgt für 5, hiermit der Wert 
worin ce den konstanten Wert 7 b; bedeutet, da b; ebenso wie e; und e,t, von un- 


abhängig ist. In Abb. 5 sind die Projektionen auf die r I-Ebene und » “-Ebene dargestellt. 
Wählt für a die Abzisse des Schnittpunktes der Geraden I und o!, so wird a 8, Q; von der 


Lage von o! unabhängig und db; nach Gl. (7) proportional dem Werte «8; 07; aus der Abb. 5 
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erkennt man leicht an Hand der zu zwei verschiedenen Werten F; und F'/’ eingetragenen 
Konstruktion, daß sich die Gerade 0) um den Punkt 0), d. i. um ihren Sehnittpunkt mit 
r; drehen muß, wenn sieh F; ändert. In dem » +1 dimensionalen Raum dreht sich also 


o' um O0! und deshalb in jeder beliebigen Projektionsebene das Bild oh um das Bild 0). Die 
Abszisse von sei v = ey. 


RA17326] / 


\ Symptore 


Ahbh. Abh. 6. 


4. Stabdehnungen nach Erreichen der Traggrenze eines Stabes. Es ist möglich, daß 
der schwächste Stab, etwa der Stab I, welcher zuerst an die Traggrenze gelangt ist, über- 
haupt keine Kraft aufnehmen kann, bevor noch ein zweiter Stab an die Traggrenze kommt, 
also bevor noch die nächste Polygonecke erreicht wird. Wir sind aber mittels unseres 
Polygones sofort in der Lage, den Sicherheitsgrad des noch übrigbleibenden statisch be- 
stimmten Fachwerkes anzugeben. Wenn das Polygon in die vi-Ebene projiziert wurde, 
braucht man bloß nachzusehen, wo die zu !; Parallele durch den Ursprung, also r; das 
Polygon schneidet, und die Abszisse dieses Schnittpunktes gibt dann den gesuchten Sicher- 
heitsgrad an. Denn wenn der Stab ! keine Kraft mehr aufnehmen kann, also 5, (v)=V 
wird, muß ®’ mit 2 (r) zusammenfallen, da ja dann mit A) auch $; verschwindet. Es fallen 
also alle mit den demnach auch s; mit zusammen. 

Es kommt also darauf an, jenen Wert » 7 zu bestimmen, bei welchem der Stab I! die 
Bruchdehnung erreicht. Die Längenänderung des Stabes /! erhält man am einfachsten aus 
der bekannten Beziehung 


le; 


worin ‚le; die Längenänderungen der einzelnen Stäbe bedeutet. Hieraus ergibt sich wegen 
so daß die Summe über alle Stäbe mit Ausnahme des Stabes ! sich erstreckt. Da voraus- 
setzungsgemäß nur der Stab 7 die Traggrenze erreicht hat, alle anderen Stäbe rein elastisches 
Verhalten besitzen, wird 
E F'; 

und daher wird «das Verhältnis y7, der Stablängenänderung an der Bruchgrenze zu jener an 
der Traggrenze 
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worin b; und a; dieselbe Bedeutung wie in Gl. (6a) besitzen und E den Elastizitätsmodul 
bedeutet. Hiermit ergibt sich für den Wert » = u,, bei welchem der Stab 1 überhaupt keine 
Kraft mehr aufnehmen kann, die Gleichung 


Sonach hat man bloß zu der Geraden 0) eine Parallele 7) zu ziehen, welche auf der Abszissen- 
achse den Abschnitt g,b; besitzt. Der Schnittpunkt von 7) mit der Geraden /, hat dann den 
gesuchten Wert u, als Abszisse. Auch diese Konstruktion kann in jeder beliebigen Projektions- 
ebene durchgeführt werden, wie dies in Abb. 4 gezeigt ist. 

Ändert sich die Fläche F; eines beliebigen anderen Stabes i, so dreht sich auch die 
Gerade 7’ um einen Punkt. Die Abszisse dieses Punktes 7’, welcher auf vr’, bzw. dessen 
Bilder Pl, auf den Bildern r), liegen, beträgt e,yr. 

Wir wollen jetzt den Zusammenhang zwischen den Werten von «7, und der ver- 
änderlichen Fläche F; des Stabes i bei unveränderlicher Fläche F, feststellen, indem auf 
den zu den Werten von F; gemäß der Gleichung T;—=F;t,; gehörigeh Geraden ti jeweils 
der Punkt mit der Abzisse u, angegeben wird. In Abb. 6 ist die Projektion auf die 
vi-Ebene dargestellt, so daß u, jeweils auf dem zugehörigen i; liegt. Zunächst sind die 
beiden zu der konstanten Fläche F7} gehörenden Geraden 1; und o; eingetragen, die nach 
früherem ermittelt worden sind. Der Punkt 0, durch den das Geradenbüschel der 0) gehen 
muß, ist hier wie folgt bestimmt worden. Zuerst ist die Gerade s; durch Verbinden des 
Punktes mit den Koordinaten vr—1 und $; , das ist der Stabkraft 8; für FH >x und 
F,>x& bei der gegebenen Belastung, mit dem Ursprung eingezeichnet. Dann ist s;, welches 
zu den Werten F;>x und F, gehört, als Parallele zu 0; erhalten. Zieht man nun dureh 


den Schnittpunkt von s; mit 1; die Parallele zu s; , so ist dies bereits jenes 0, das zu 
F;>» gehört. Denn auf diesem ot müssen sich bei konstantem F;>x jeweils bei 
veränderlichem F7, einander zugeordnete Werte von s; und 1; schneiden, und wir haben 
dureh zwei Paare solcher Werte, nämlich s; und 1; für das gegebene F), und s; für 
F}A>x die Gerade 0} festgelegt. Ihr Sehnittpunkt mit r', welches hier mit der »-Achse 
zusammenfällt, ist der gesuchte Punkt durch den alle für verschiedene F; 
eehen müssen. Durch den Punkt pP! mit der Abszisse y,ec; geht dann das Geradenbüschel 


der zu 0} parallelen 7. Man kann sich aber das Zeichnen der beiden Geradenbüschel der 
0! und r ersparen, wenn man durch den Schnittpunkt von 1 und /; die Parallele & zu s; 
zieht und durch deren Schnittpunkt N: mit r ein Geradenbüschel £; parallel zu den zuge- 
hörigen s; legt. Die Schnittpunkte von 7; mit £; sind dann mit denen der r identisch; pro- 
jiziert man sie schließlich auf die zugehörigen i;, so erhält man eine gleichseitige Hyperbel, 
deren Koordinaten einander zugeordnete Werte von a; und T;, also auch F; sind. Die Assym- 
ptoten der Hyperbel können, wie aus Abb. 6 ersichtlich, leicht gefunden werden. Der 
Schnittpunkt von I; mit r} ist ein Punkt der Hyperbel. Da der Stab I nur dann als erster 
Stab an die Traggrenze kommt, wenn 7 über dem Sehnittpunkt von I; mit 0 liegt, so kommt 
nur jener Teil der Hyperbel in Frage, bei dem die Ordinaten der Hyperbel größer sind, als 


die Ordinate T;, des erwähnten Punktes. 


5. Die Abhängigkeit des Sicherheitsgrades von den Stabquerschnitten. Wir sind nun 
in der Lage, den Sicherheitsgrad eines einfach statisch unbestimmten Fachwerkes als 
Funktion der Fläche eines Stabes anzugeben. Es kommt häufig vor, daß ein solches System 
durch Weglassen eines Stabes als statisch bestimmtes Fachwerk berechnet wird. Es ist 
nun die Frage, ob durch das nachträgliche Einfügen eines Stabes der Sicherheitsgrad ver- 
größert oder verkleinert wird. 

Wir wollen annehmen, daß die Fläche F, des Stabes p von Null bis Unendlich wächst und 
verfahren wie folgt, um den Sicherheitsgrad zu bestimmen, welcher zu einem beliebigen 
Werte von F,, gehört. Zunächst zeichnen wir die Projektion des Polygones in die v p-Ebene. 
Trägt man die Gerade 0), ein, so kann man nach früherem für jedes F, die dazugehörige 
Lage von s, angeben, da ja F, der Tragkraft T, bzw. T,’ proportional ist, und man erhält 
die noch fehlenden Polygonseiten p, und p,', als Parallele zur »- Achse in den Abständen 
T,-- F,t, und T’=F,t,. Man erkennt, daß der Stab p dann als erster Stab an die Trag- 
grenze kommt, wenn die Gerade 06’ im Inneren des Polygones liegt. Die Koordinaten der 


Polvgonseite m, werden daher in dem in Abb. 7 dargestellten Fall einander zugeordnete 
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Werte von T, bzw. F,, und des Sicherheitsgrades darstellen. Allerdings ist es möglich, daß 
sich der Stab schon früher ausschaltet, wie dies in Abb. 7 zwischen den Werten F,,=0 
und F,, der Fall ist, da hier die Gerade r/, ebenfalls im Inneren des Polygones liegt. Dann 
kommt der Sicherheitsgrad des übrigbleibenden statisch bestimmten Systemes, welchem der 


Stab p fehlt, in Frage, der also bis zum Schnittpunkt von 7), mit m, konstant bleibt. Erst 


für Flächen zwischen F,, und F,,, ist der Sicherheitsgrad durch die Abszissen der Polygon- 
seite m, dargestellt. In diesem Bereiche tritt Unbrauchbarkeit des Systemes in der Weise 
ein, daß zuerst der Stab p. dann der Stab m an die Traggrenze gelangt. 

Wächst die Fläche F, weiter, so kommt als erster Stab der Stab m an die Traggrenze, 
als zweiter zunächst, d. i. für Werte von F\, zwischen F',,, und F,,, der Stab p; dabei wechselt 
die Stabkraft S, (v») jedoch ihr Vorzeichen. Der Sicherheitsgrad des Systemes ist jetzt durch 
den Schnittpunkt von p,’ mit m, gegeben, welcher aber aus Zweckmäßigkeitsgründen auf die 
Gerade p,. die zu demselben F, wie p,' gehört, projiziert ist. Dies gilt bis F,,; von da an 
bleibt die Sicherheit bis F,, konstant, da der Schnittpunkt von m, mit n, nunmehr für den 
Sicherheitsgrad maßgebend ist. Das System wird hier dadurch unbrauchbar, daß zuerst der 
Stab m, dann der Stab n an die Traggrenze gelangt. Wächst schließlich F, noch weiter bis 
zu einem unendlich großen Wert, hat also schließlich s, die Riehtung von 0), so stellt die 
Abszisse des Schnittpunktes von 4, mit m, den Sicherheitsgrad vor, da jetzt zuerst der Stab 
y, dann der Stab m die Traggrenze erreicht. 


i 
| 
| 
| 
| 
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Abb, 7. 


Es ist noch möglich, daß der Stab m und der Stab y an die Bruchgrenze kommen, be- 
vor noch der zweite Stab die Traggrenze erreicht hat. Um sich bezüglich des Stabes m, bei 
welchem dies in dem Intervall zwischen F',, und F',, der Fall sein kann, hierüber Klarheit 
zu verschaffen, ist die früher erwähnte Hyperbel für diesen Stab gemäß dem in Abb. 6 dar- 
gestellten Verfahren einzuzeichnen. Man sieht, daß in dem vorliegenden Falle durch das Er- 
reichen der Bruchdehnung des Stabes m in dem ganzen Intervall von F,,, bis F,, der Sicher- 
heitsgrad heruntergesetzt wird, da «das statisch bestimmte Fachwerk, das durch Ausschalten 
des Stabes m entstanden ist, eine kleinere Sicherheit besitzt, die durch die Abszisse des 
Punktes 9 gegeben ist. Es ist aber möglich, die Belastung auf jenes Vielfache » zu steigern, 
das durch die Abszissen der Hyperbel gegeben ist, und die Koordinaten der Hyperbel stellen 
also einander zugeordnete Werte des Sicherheitsgrades » und der Fläche F, vor. Eine der- 
artige Verminderung des Sicherheitsgrades infolge der Bruchdehnung kommt im allgemeinen 
nur bei Druckstäben in Frage, da für Zugstäbe der Wert von g viel zu groß ist, als daß die 
Bruchgrenze des ersten Stabes vor der Traggrenze des zweiten Stabes erreicht wird. Dies 
ist auch der Grund, warum statisch unbestimmte Systeme in den allerseltensten Fällen durch 
Zerreißen eines Zugstabes zugrunde gehen. Man sieht aber auch weiter aus dem in Abb. 7 
dargestellten Beispiel, daß die Vergrößerung des (uerschnittes eines Stabes nicht immer eine 
Erhöhung des Sicherheitsgrades zur Folge haben muß. Im unserem Falle erreicht man den 
größten Sicherheitsgrad dann, wenn man den Stab p überhaupt fortläßt und es ist daher 
der oftmals geübte Vorgang, ein statisch bestimmtes System durch Einfügen eines weiteren 
Stabes verstärken zu wollen, nicht immer von Erfolg begleitet. 173 
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Der ebene Spannungszustand. 
Von Karl Hager in München. 


m Jahrgange 1930 dieser Zeitschrift hat Dr. Craemer in einem Aufsatz „Spannungen in 
| heran Balken bei feldweise wechselnder Belastung* eine Spannungsfunktion für den 
ebenen Spannungszustand gewählt und damit die Unzulänglichkeit der Navierschen Biegungs- 
theorie für hohe Tragwände nachgewiesen. Da nun diese Funktion für die Rechnung nicht 
sehr bequem ist, soll in folgendem der Versuch gemacht werden, eine einfachere Form für 
die Biegungsspannungen zu finden, welche aber doch die hinreichende Genauigkeit für den 
praktischen Gebrauch noch gewährleistet. 

Wenn & und £ die Verschiebungen eines Punktes (, 2) sind unter der Einwirkung der 
Spannungen 9%, T, so bestehen bekanntlich folgende fünf Gleichungen (Abb. 1). 


, Ö | 
00, , OT;x_ Q 
0; £ Ox (2) 
05,0 
Aus (1) kann man, wie angedeutet, 7,- eliminieren. 
N? 06, 
.x° 02 
Or 
Aus (3) erhält man +3 und aus (2) daß sich ergibt 
OT l 1 00, 1 0°, E 
0% 
| 
| 
Tız | 
| 
n| | 
R 
zz 
Abb. ] Abb. 2 
m 
;, 
aus (2 = setz 
‚(2m+l 5 
Mit Hilfe von (I) erhält man dann die zweite Hauptgleichung 
2.-—— | (Il) 


Aus (I) und (II) ergeben sich die Differentialgleichungen für die beiden Normalspannungen zu 


0’ 0. o. 

- +2. | 


202 
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0 
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ls sind die Differentialgleichungen der biharmonischen Funktionen, welche auch bei der 
Berechnung ebener Platten benutzt werden'!). 
Die Gl. (III) wird erfüllt, wenn man setzt 


0, =Ü-.cos 


] ] z+2-(, 


wobei der Koordinatenanfang in die Mitte des unteren Balkenrandes gelegt ist (Abb. 2). 
Für freie Auflagerung ist in den Punkten # + ,„ 0.=0. Auch diese Bedingung wird 
von dem angenommenen Ausdruck erfüllt. 


0% 
( 


Mit Hilfe von (I) führt diese Gleichung also auch zu einem Ausdruck für die lotrechte 
Spannung ©.. 

Jedoch soll zunächst o, weiter betrachtet werden. 

Außer der Differentialgleichung (III) muß der Wert o. noch zwei weitere Bedingungen 
erfüllen, welche sieh ergeben, wenn man den Balken lotreceht abschneidet und die inneren 
Kräfte zur Herstellung des Gleichgewichtes wieder ansetzt. Unter der Voraussetzung nur 
lotreehter Lasten müssen somit die beiden Gleiehgewiehtsbedingungen an diesem Querschnitt 
erfüllt sein, welche sieh in den Gleichungen ausdrücken 


h h 

\o,.dz2=0 und 


Wenn die Gleiehung für alle Werte . gelten soll, müssen die Faktoren von x Null 


en 
werden und ebenso das konstante Glied. Für cos kann man setzen |1 — | ] 
Die Konstante der Gleichung ist 


Man erkennt, daß die beiden Gleichungen nicht übereinstimmen. In Gl. (4) kommt der 


Faktor , vor, und diesem entspricht in der Gl. (5) der Faktor % Es wird später auf diese 


Unstimmigkeit zurückzukonmmen sein. Vorerst soll die zweite Gleiehgewiehtsbedingung, die 
Momentengleichung, behandelt werden. 


2 
—a”)\zde 


h 

u 

h 

7 .2d2=-, Sn =H 
I I 


en 
M,=C:cos > (4 


') Vgl. Dr. Heneky, „Über den Spannungszustand in reehteekigen, ebenen Platten“. München 1913. 
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Auch in dieser Gleichung müssen die Faktoren von «°” und die Konstante Null werden, 
wenn sie für alle Werte von x Geltung haben soll. Es ist daher nötig, für die weitere prak- 
tische Rechnung einen Belastungsfall anzunehmen, der aber mit Rücksicht auf den gewählten 
Spannungsausdruck zur Balkenmitte symmetrisch sein soll. Wir wählen die gleichförmig 
verteilte Belastung p auf die Längeneinheit, und die Balkenbreite soll eins sein. 


2 
2 2? h? " 
Der Faktor von .r? ıst 
Die Konstante der Gleichung ist 


Auch die Gl. (6) und (7) zeigen eine Unstimmigkeit. Dem Faktor 7° 9869 in (6) ent- 
spricht der Faktor S in (7). Diese Unstimmigkeiten in den Gl. (b), (5) bzw. (6), (7) ließen 
sich vermeiden, wenn man den Ausdruck für o, durch eine Reihe ausdrücken würde, aber 
diese Umständlichkeit soll ja vermieden werden. Man hat daher zunächst nach dem Zweck 
der Rechnung die Konstanten zu wählen. Im Koordinatenanfangspunkt werden die Glieder 
mit .-” an sich Null. Man wird also bei Rechnungen mit den Maximalmomenten symmetrischer 
Belastungen die Gl. (5) und (7) für die Bereehnung der Konstanten € und Z benutzen, während 
man für andere Punkte des Balkens in den Gleichungen am besten einen Korrektionsfaktor 
anbringen muß. Wobei übrigens bemerkt sei, daß für die meisten Aufgaben die Unterschiede 
der 4 Gleichungen von € und Z praktisch fast bedeutungslos sind. 

Wir haben für cos 7 die beiden ersten Glieder seiner Reihe gesetzt. Man kann nun 


in dieser Annäherung den Korrektionsfaktor anbringen. 


Will man die Unstimmigkeiten ganz beseitigen, so würde in den Gl. (H und (5) y 1,27 

und in den Gl]. (6) und (7) y 1,23 sein müssen. Um nun nieht zu weit von dem Wert 

08 7 abzuweichen, wird man y 1,1 bis 1,2 zu wählen haben, wenn die Rechnung Punkte 

ın der Nähe der Auflager betrifft. 


Wir haben oben bereits gesehen, daß mit Hilfe der Gl. (I) auch ein Ausdruck für o. 
gefunden werden kann. 


0’ 0’ 0, 
02° l I 


I 


Wenn der untere Balkenrand (z - 0) unbelastet ist, muß dort 6. 0 werden, d.h. 


cos" 
Für 2. h, oberer Balkenrand. welcher die Belastung p trägt, ist 0. p 
‚_»,€ 


| 
| 
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Zur Berechnung der Schubspannung 7 muß man auf die Gl. (1) zurückgreifen. 


00, Or 
- sin 22 =— 


Da für 2-0 aueh T=0 wird, ist die Konstante T=0., 


Bei symmetrischer Belastung ist auch für #0 7-0. Auch diese Bedingung erfüllt 
die Gleichung für r. 


Es sind nunmehr noch die Konstanten € und Z aus den oben entwickelten Gl. (4), (5) 
(6) und (7) zu bereehnen. Zunächst sollen Punkte des Trägers behandelt werden, welche in 
der Nähe des Koordinatenanfangs liegen, so daß in den Spannungsgleichungen die Glieder 
mit =? an sich Null werden und nur noch die Konstanten Null werden müssen. Es sind daher 
die beiden Gl. (5) und (7) zu benutzen. 


Sie ergeben: 
p, Sin 
< 
| 


Werden andere Punkte des Balkens betrachtet, so benützt man, wie oben bereits gezeigt, 
die Gl. (4) und (6) unter Verwendung eines Korrektionsfaktors y'. Sie lauten dann: 


H-n p 

1 h 


Hieraus ergeben sich die beiden Konstanten zu 


7 


Man erkennt, daß die Werte der Konstanten Z in beiden Fällen gleich sind. 


Man kann auch für das Verhältnis der Balkenhöhe h zur Stützweite po in die 


Gleichungen einführen und erhält dann folgende Ausdrücke: 
Für Punkte in der Nähe des Koordinatenanfangs 


S-+tarZinan -Skojar’ Sinar ar) 


und für andere Punkte die mit dem Korrektionsfaktor versehenen Werte 


pr-Zinar 


/ur Berechnung der Schubspannungen muß man zu den Gl. (1) zurückkehren. 


sin , 


pr N a 


. LIT T r m 


| : | 
B: 
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Für 20 muß r=0 werden, so daß auch T=0 zu setzen ist. Bei symmetrischer Be- 
lastung zur Balkenmitte wird r—=0 für a=0. Auch diese Bedingung wird von der Gl. (8) 


erfüllt. Dagegen wird die Bedingung nicht erfüllt, daß 7 auch 


für den oberen, belasteten Balkenrand Null wird. Vielleicht ließe 
sich das in einem komplizierteren Ausdruck erreichen. Aber auch 
in der eingangs erwähnten Arbeit, welche als Spannungsfunktion 
einen sehr komplizierten Ausdruck in Gestalt einer unendlichen 


Reihe benützt, ist die Bedingung nicht erfüllt. 


Zum Schlusse seien die praktischen Ergebnisse der vor- 
stehenden Entwicklungen noch an einem Zahlenbeispiel gezeigt. 


Zahlenbeispiel: 


a=10; 100 t/m?’= 10 kg/em?. 
Die Last soll oben liegen (Abb. 3). 


Nach (5) und (7) Abb. 3. 
C=-163; Z=—-%2 
— iD 
438; z=h 0.=+12369 Navier 0,=+150; 126.9 
h 30 
o.—=5150 Navier 00 
2 515 


Nach (6a) und (7a) 


C=151; Z=—236,2 
63 
| h 3 | 
-021; s=— 0,=389;, h 0.=66,8 
z=04l; z2=0, = —-16,35; =h; 0,=+33,8 Navier 27; 
h > 
c=04l; 2=— 0,=93W0; 2= 4 h o0.=43! 
0,51: und z—h Öx = Navie 0 
h 
s=— r=Bl2 Navier 
Dasselbe Beispiel mit abnehmendem Höhenverhältnis «a 
| 
2180;  Z= - 1190 
_ 672 
z2=0 0,=—-500; z=h 0,=+510 Navier 672; s10 0,83 
| | 
h=030 a= 0 =253720; Z= 1146535 
1460 
2=0 0%= z=h 0,8375  Navier 460; 0,54. 


0,59 


Man erkennt, daß mit abnehmendem Höhenverhältnis die Randspannungen sich dem 


Navierschen Werte nähern. 


In der Abb. 3 sind die wichtigsten Rechnungsergebnisse für eine quadratische, gleich- 


förmig belastete Wand als Träger auf zwei Stützen zusammengestellt. 


Würde die Last an dem unteren Balkenrande angreifen, so wäre nur der Koordinaten- 


anfang in die obere Balkenkante zu legen. 


Zusammenfassung. Ausgehend von den fünf Differentialgleichungen, welche für den 
zweiachsigen Spannungszustand zwischen den Spannungen und Verformungen in jedem Körper- 
punkte bestehen, werden für die Spannungen o,, o. und r Gleichungen entwickelt, welche 


die Differentialgleichungen und die Randbedingungen erfüllen. 
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Zur Statik eingespannter Rechteckplatten. 


Von Harry Schmidt in Köthen. 


De Benutzung eines früher!) angegebenen Integrationsverfahrens habe ich kürzlich?) die 
Durehbiegung einer frei aufliegenden Rechteckplatte bei statischer, gleichförmiger Be- 
lastung eines parallel zu den Plattenrändern abgegrenzten Plattenteils berechnet, womit die 
Mehrzahl der praktisch wichtigen Lastverteilungen durch geeignete Spezialisierungen erfaßt 
werden kann. In der vorliegenden Mitteilung soll gezeigt werden, daß sich der gleiche Be- 
lastungszustand auf genau entsprechende Weise zunächst für die zweiseitig eingespannte, 
zweiseitig frei aufliegende Rechteekplatte behandeln und alsdann der Fall der vierseitig ein- 
gespannten Platte sich auf das Problem der Auflösung von zwei Systemen je einfach unendlich 
vieler Gleichungen mit einfach unendlich vielen Unbekannten reduzieren läßt. 


1. Bezeichnungen und Ansatz. Die Plattenmittellläche möge das durch und 
0°. 9°. 5b gekennzeichnete Gebiet B der X Y-Ebene bedecken; für bzw. a sei feste Ein- 
spannung, für y==0 bzw. b hingegen freie Auflagerung angenommen, und die zu betrachtende, 
auf die Flächeneinheit bezogene Lastverteilung p (#, y) werde mit festem p, durch den Ansatz 
in B,. 


beschrieben, wobei R, das aus Abb. I ersichtliche rechteckig berandete Teilgebiet von B be- 
deutet. Setzen wir sodann 


»(2,y)= (1,1) 


m’ Eh’ 


12(m? — 1)’ 
1Y 
d 
| Bo 
Abh. 1. 


unter h die Plattendieke sowie unter E bzw. m den Youngschen Elastizitätsmodul bzw. die 
Poissonsche Querkontraktionszahl des Plattenmaterials verstanden, so ist die Durchbiegung 
(we, y) als dasjenige Integral der linearen inhomogenen Differentialgleiehung 


l 


zu ermitteln, das den acht Randbedingungen 


umel 
(el, — (ir, \ 


genügt. Durch Heranziehung der bekannten, für jedes positive «a, gültigen Relation 


en für? 0% (1,5) 


läßt sich aber (1.1) durch 


I!) HM. Schmidt, Zs. f. ange. Math. u. Mech. ®, 401, 1929. 
Schmidt, Zs. f. Phys. 68, 423, 1931. Vgl. überdies die inzwischen erschienenen Arbeiten: Ingenieur-Archiv 
>. 449, 1931: Zs. F. ang. Math. u. Mech. tt, 426, 1931: Zs. f. Phys. 74, 130, 1932. 
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ie mit 
1e Po ((e‘“ xı)2 — Xp)? 
. . . . 
. ersetzen, so daß wir unter Berufung auf das lineare Superpositionsprinzip für die gesuchte 
h Durehbiegung w(#,1y) zunächst den Ausdruck 
Po r - 
dl 
mit 
erhalten, falls 
gesetzt und mit (#,9y; 2,{) diejenige Lösung der Differentialgleichung 
bezeichnet wird, die den (1,21) und (1,22) entsprechenden Randbedingungen 
und 
[u (se, [u (#, 9; 2, (1,6 
genügt. 
1 


2. Lösung des Randwertproblems. Um nun das Randwertproblem (1,6) mit (1,61) und (1,62) 
zu lösen, gehen wir von dem partikulären Integral 


der inhomogenen Gleichung (1,6) aus und bestimmen eine der homogenen Gleichung 


genügende Funktion von der Form 


was 
2.) = Alz,{) sin (yz2)—+ biz,{) cos(y2) 


+ y- cos (yz) 


liefert. Die hierin auftretenden Koeffizienten wählen wir derart, daß die Summe 


u, (2,952, )+ uw (a, + f(y; 2,{), 


die Randbedingungen (1,62) erfüllt, wodurch sich 


[2 
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22. sin?’ (bz 
2’ 
Ba, C(2,{)= 2,22) 
Die -(2?+L°)- cos 
12,6) 2z-sin (bz) 
ergibt. Setzen wir nun die Sinusreihe 
I K,(e,d)- sin (2,3) 
mit 
v7 9% 
Pı % (2,31) 
an, so wird | 
h 
2 | 
oder 
2 ß, 1 — (- 1)”. eb» 
und da die Funktion 
ur K,(2,{): (x; sın (P, Y) (2,4) 
mit 
(#52) = c,” (2) Sin (e — a)] + (2) - Sin 
| (2,41) 


+ c,® (2)- x - Sin [P, (& 


der homogenen Gl. (2,11) sowie den Randbedingungen (1,62) genügt, so muß der Ansatz 


(2,9; 2,0) 
oder 


2,5) 


e,® (2)- (x — a): Sin 


SL 
jet? + sin (P, 9). 


2,51) 


das sämtlichen verlangten Nebenbedingungen entsprechende Integral von (1,6) darstellen, 
falls die Koeffizienten von (2,41) so bestimmt werden, daß 


eilt. Dies führt aber zu den Gleichungen 


— 


z2+P, (Pf, a), + Pr (2) 


die dureh 
| 


Zin(P,al 


elz.19,a2- zin 4) 


Pr 


(2) 
Sin? (ß,a) — 
Pr zın (Pr dt) 
(2) 


(pP, a) 


zelöst werden. 


(2) 
Zin (pP, a) 
Zin(p,a)- a) Sin’ (P,a)+ 


(P,a)-c,"” (2) 


Zin(P,a)' 


a, 
Zin (pP, a) 


Pi 


2. Sin? (P, 
a) 
- Zin (pP, a) 


pP, az- Zin(P,«a) 


= 


2,6) 


Ay 
| 
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3. Auswertung der komplexen Integrale. Zur Auswertung der Integrale (1,52) und (1,51) 
benutzen wir eine früher”) bewiesene Verallgemeinerung von (1,3), die durch 


2ni 
Ay 
0 für 
mit 
So 
A,” 
F iz) — (2 — 2, x o) (6 


zum Ausdruck gelangt, wobei k(z\ für z>w als beschränkt vorauszusetzen ist, und wobei 


sämtliche (isoliert liegenden) Nullstellen z= 2, der analytischen Funktion Fiz) Realteile kleiner 
als a, 0 besitzen sollen. Wegen der in Analogie zu (3,11) bestehenden Entwicklungen 


| I 4 
| 
22T 
sowie 
3 
l 
l 
22” 
erhalten wir hiernach gemäß (2,1) 
| er für y> | 
% — (3,2) 
2 ni 
| für y<y, 
mit 
(320 


‚cos[y 


während sich gemäß (2,2) nach Einführung von (2,21) und (2,22 


mit 
sin(y2) boy | 


- 


b-cos(bz)-sin(yz2) y-sin(bz)-cos(yz) f 
A 


| 
2z-sın’ (bz) 


cos |(b | 


s), H. Schmidt, Zs. f. Phys. 39, 474, 1926 
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sowie gemäß (2,4) und (2,32) 


+ im 
mit 
(-1P-sin 
F.(x,9; 2) Byb (2 - Be (3,41) 


berechnet. Durch Eintragung von (3,2), (3,5) und (3,4) in (2,5) ergibt sich sodann auf Grund 
von (1,52) 
mit 
oder wegen (3,41) 


\ sın (P, 4 2) 


sowie mit 
für O<y<sy, 
((y: 2) G,(y;2) für 
G,(y:2) für | 
wobei 
9.19: 2), | 
G,(y:2) = - 9.(y; 2), 


zu setzen Ist. 
4. Formulierung des Endresultats. Da gemäß (3,21) und (3,31) 


sin (bz) sin (42) — cos(y,2)- sin |(b )z] 
sin (bz) 


(9,2). cos|(b -cos(yz)ı sin(y,2)- sin|(b-)z] 
22". sin (bz) 
(y,„2) eos (bz)} - sin (y2) 


22°. sın? (bz) 


+ 


gilt, so sind nach Eintragung von (3,5) in (1,51) die Voraussetzungen zu erneuter Anwendung 
von (3,1) sämtlich erfüllt. Zunächst ergibt sich wegen (3,52), (3553), (3,2D und (3,31) 


falls 
mt 
sowie 
(y) für (0): Y, | 
für yo | 
mit 


1 


rw 


7 
- 
( 


1) 


or 
> 
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gesetzt wird. Benutzt man an Stelle der geschlossenen Ausdrücke (2,1) und (2,2) die Reihen- 
darstellung (2,3), so wird man durch Eintragung von (2,5) in (1,52) und Anwendung auf (1,51) 
zu dem Ausdruck 


[9] 


geführt, der nichts anderes als die Fourierentwicklung von (4,3) darstellt, und der unmittelbar 
die Berechtigung der in (4,1) hinzugefügten Gleichheitszeichen erkennen läßt. Ferner wird 
gemäß (3,51), (2,41) und (2,6) 


Ag | 
Aan—ix 
wobei 
v—| 
mit 
4,5 
— ax-Sin(P,xr)- Sin(P,a) 
und 
_ 
l 
sowie mit 
= za P,-Sin[p,(2 —a)] —a)- . . (4,53) 
und 
1 | 
2 
(4,54) 
- Sin[?,(a— ,)]} 


| 
m 


y) 
7) 
X; q 
Abb. 


gilt. Die Einführung von (3,5) in (1,51) liefert demnach gemäß (4,1) und (4,4) das in Abb. 2 
veranschaulichte Endergebnis 


W,(e,y) für 
We, Wr, y) für ra, 
das sich auch in Form der Doppelsinusreihe 
Wie, y) —) Las 2 
mit 
a 


+: 
ch. 
nd 
>) 
2) 
| 
10* 


Ztsehr. f. angew. 
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und 
a h 
4 
schreiben läßt. Diese Koeffizienten Lu, berechnen sich am einfachsten dadurch, daß man 
(2,51) nach Vornahme der Fourierentwieklung der Funktion e&®—+9,(“:z) in (1,52) einträgt 
und die hiermit gebildete Relation (1,51) gemäß (3,1) auswertet; auf diese Weise ergibt sich 
die auch direkt zu bestätigende Beziehung | 


“> 


Au Pr" . 4 Pr) 


falls 
kur? — Sin 1)" 


gesetzt wırd. 


5. Spezialfälle. Wird das Belastungsgebiet B, (vgl. Abb. 1) mit der vollen Platten- 
mittelfläche identisch, so ist sowie und b zu setzen, und (4,6) ver- 
einfacht sieh zu 


4 
(3,1) 
| 
mit und (1,59) sowie mit 
Sin(p,a)y, 
während in der Darstellungsform (4.7) 
A (au . . . . . . . . . . Im 
mit 
| (- 1)" -Ein(p,a) + P,a 521) 
wird. 


20 
übergang 90, so ergibt sieh für die Plattendurehbiegung unter dem Einfluß einer gleich- 
förmigen Linienbelastung vom Betrage P, pro Längeneinheit, die längs der zu den beiden 
Plattenrändern und a parallel gelegenen Strecke y, angreift, gemäß 
(1,5) und (4,6) der Auselruck 


Wählen wir ferner, Sowie und vollziehen den Grenz- 


0 2 


| W,(x,9) für | 
wobei 


v=] 


= 
» . 
| 


1939 Schmidt, Zur Statik eingespannter Rechteekplatten 149 


mit 


gilt, und wobei W,(r,y) gemäß (4,5) mit (4,5D und (4.53) sowie mit 


Pr - 
und 


zu bilden ist; die Fourierentwieklung (4,7) ergibt 


ab By: (au: + (5,9) 
nit 
| (— 1" - Sin (f,a) — a 
Um schließlich auch den Fall einer im Punkte (#,., »,) angreifenden Einzellast vom 
Betrage P, zu erfassen, haben wir in (5,3) und (5,31) die Substitution P, = zu machen, 


9 sowie 9,9, +9 zu setzen und den Grenzübergang 9>0 zu vollziehen, was für 
die Durchbiegung 


P, 


mit 
liefert, falls 
mit 
sowie 
W, Sin? — sin (P, 9) 
v=] 


mit (4,51), (4,53) und (5,42), (9,45) gesetzt wird; die Fourierentwieklung (4,7) erhält alsdann 
die Koeffizienten 


sın (P, Tun 


(2.9) 
mit (5,51). 


6. Die allseitig eingespannte Platte. Uns nunmehr dem Fall der vierseitig eingespannten 
Rechteckplatte zuwendend, gehen wir von den Funktionen 


fa” y- b)]- sin (au) sin sin (PM). 
v—] 


(6,1) 


fu y) = (y b)- Sin (au - sin (aux) sin (au): - sin 


ech. 
2) 
an 
igt 
ch 
1) 
) 
T- 
| 
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mit | 
Pr 


(6,10) 


sowie 


ui 
(62 
F, (x a)- Zin(p,x)- sin (P,y) = sin N Bu® - sin (au) 
u J 
mit 
taup, 
Au a: (a + (rar, 
(6,21) 


aus, bilden den der homogenen Gleichung 
(ar, y) 


eenügenden Ansatz 


SF. 


= I (a, y) + fu? (a, 
— ] 


+ 310,” + F,”(&,y)} 


und suchen hierin über die willkürlich wählbaren Konstanten derart zu verfügen, daß neben 
den bereits erfüllten Beziehungen 

| = 0 — | W," (r, | W,* (r, Yu =0 — | 0 (6,31) 
überdies die beiden Randbedingungen 


=0 


befriedigt werden. Indem wir (6,3) mit Benutzung von (6,1) in der Form 


W*ay)= 2 | Sin[ß,(& a)]+C,® --a)- Sin (B,x) 
1 


schreiben, liefert dies die beiden Relationen 


| 


aus denen sich 


\ 
Y ) 
(P, a) + (— 1)" - Pr a, Au Cu + Bus‘ Cu”) 
(' 2 
Pr (6,38) 
( 1)" - Au om (Pr 1)" a, Au Cu” + Bus 
m 1 
Sin? (P, a) — a? 


| 
| 


lech. 


1m) 


1) 


en 


{ 3 H ft r N 
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hereehnet. Greifen wir sodann auf die durch (1,5) eingeführte, in (4,6) bzw. (4,7) angegebene 
Funktion W (x, 9) zurück, so wird der Ausdruck 


die Durcehbiegung der allseitig eingespannten Rechteckplatte darstellen, falls die beiden in 
(6,3) noch frei verfügbaren Konstantensysteme (C.“” und C,” derart bestimmt werden, daß 


104 yab 
gilt. Dies gelingt aber durch den gemäß (6,3) und (6,2) gebildeten Ansatz 


W*(, Zinl[au(y— b)| + Cu” (y b)- Sin (au y) 


demzufolge mit Benutzung von (4,7) 


) | 


Cu zin (Y b) + du Co] (Y 21 
Cu” zin (An Au (Y h) Co] (au Y) 


wird, so daß (6,41) die Gleichungssysteme 


> Pr zın (au h) Cu’ — aub- Cu” 


— 


Pr‘ (0,9. Au? + 
(' 
(-1)"- taub: Cu” + Sınlaub)- Cu” 


+ + 0,8. Buy 


ergibt. Hieraus folgt nach Einführung von (6,33) unter Beachtung von (6,11) und (6,21) 
durch Addition 


. 1 + 1)’ Cu . (an bh) Au h} 


— v 
es) 
Sin?’ (P,a) a? 
1 


=0, 
hingegen durch Subtraktion 


41 ( 1)", Cu?) (au h) + Au 


1 


Sin’ (d,a) Pre 


o=1 


womit zwei voneinander unabhängige Systeme von Bestimmungsgleichungen zunächst für 
Cu” bzw. Cu” + Cu” und damit auch für und Cu“ aufgefunden sind. 175 


| ' 
) 
| 
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ur | 
|) | 
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Die Berechnung der Gleichgewichtslage von gemessenen 
Schwingungen auf Grund der Fehlertheorie’. 


Von Max Schuler in Göttingen. 


enn ich eine so einfache Sache wie die Berechnung der Gleichgewichtslage von 
N Schwingungen hier behandle, so geschieht dies deshalb, weil man in allen Büchern 
falsche Darstellungen darüber findet. Ich glaube, daß ein kurzer Hinweis genügen wird, um 
die ererbten Fehler in Zukunft auszumerzen. 

Kohlrausch sagt in seinem Lehrbuch der praktischen Physik über die Bestimmung 
der Gleiechgewichtslage: „Man wähle eine ungerade Zahl von Umkehrpunkten, nehme das 
arithmetische Mittel aus den Beobachtungen auf der einen Seite und aus denen auf der 
anderen Seite, Aus diesen beiden Zahlen bilde man wieder das Mittel.“ Diese Vorschrift, 
die man fast in allen Lehrbüchern findet, ist falsch aus folgenden Gründen: 


I. Es ist unzweckmäßig, bei einer geraden Zahl von Beobachtungen eine Beobachtung 
herauszuwerfen, die man vielleicht mit vieler Mühe gemacht hat. 

> Die eine Seite wird mehr bewertet als die andere Seite; denn das Mittel der beiden 
Seiten ist aus einer verschiedenen Anzahl von Beobachtungen gewonnen. 

3. Kohlrausehs Rechnung gibt falsche Werte bei gedämpften Schwingungen. 

I. Man erhält trotz vieler Beobachtungen keinen mittleren Fehler und damit keinen 
Aufschluß über die Genauigkeit der bestimmten Gleichgewichtslage. 


Im Gegensatz zu Kohlrausch sagen die Herren vom Vermessungswesen’): Man be- 
stimme diejenige e-Funktion, die sich den beobachteten Umkehrpunkten am besten anpaßt, 
für die also gilt: 

> — Minimum. 
Diese Vorschrift macht nieht nur viel Rechenarbeit, sondern ist ebenfalls falsch, wie ich im 
folgenden beweisen werde. 


Bei der Berechnung der Gleichgewichtslage wollen wir drei Fälle unterscheiden: 


a) Keine Dämpfung: d. h., während einer einzigen Periode kann der Schwingungs- 
abfall vernachlässigt werden. 

b) Kleine Dämpfung: d. h., während einer einzigen Periode kann der Schwingungs- 
abfall als geradlinig angesehen werden. 

ec) Große Dämpfung: d. h., man muß der Rechnung eine e-Funktion zugrunde legen, 
welehe die Amplituden begrenzt. 


Haben wir » Umkehrpunkte beobachtet, so ergibt dies an » 1 Schwingungsbogen. Im 
Falle daß keine Dämpfung vorhanden ist, bilden wir das Mittel /, jedes einzelnen Schwin- 
gungsbogens (siehe Abb. 1): 


da, 
Als Gleichgewiechtslage ergibt sich: 
a n 
Aa 
= — 
y . . . . . . . . . . . . . . 


Abb. 1. Falla: keine Dämpfung. Abb. 2. Fall b; kleine Dämpfung. 


I!) Nach dem bei der Tagung des Aussehusses für mechanische Schwingungen des V.D.I. in München 1931 
gehaltenen Vortrag „Berechnung der Gleichgewiehtslage von Schwingungen“. 
:), Basch, Mitteilungen aus dem Markscheidewesen 1917 und P. Wilski, ebenda 1922. 


\ 
| 
dl 
| 
| 
7; (L > 
A; 
Zu | 


h. 


n 
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Wir haben das arithmetische Mittel aus sämtlichen Umkehrpunkten zu 
nilden, wobei aber der erste und letzte Umkehrpunkt nur halb zu be- 
werten sind. Diese Bewertung ist auch logisch berechtigt, da man Störungen in der Nähe 
es ersten und letzten Umkehrpunktes nicht so gut erkennen kann wie Störungen innerhalb 
des Sehwingungsbildes. Bei dieser Rechenregel ist es gleichgültig, ob die Zahl der Umkehr- 
punkte gerade oder ungerade ist. Man hat stets gleiche Gewichte für jede Schwingungsseite. 
Die Rechenregel ist jetzt einwandfrei und auch einfacher als die von Kohlrausch ange- 
vebene. Aus den Abweichungen » der Bogenmitten L gegen die Gleichgewichtslage @,„ 
erhält man den mittleren Fehler für eine Bogenmitte L: 


'2v / 2v 
=rt + .) 


und den Fehler der Gleichgewichtslage: 


Dieser mittlere Fehler ist riehtig, wenn jeder Schwingungsbogen eine „unabhängige Be- 
obachtung“ ist. Dynamisch ist dies sicher der Fall. Liegen aber Ablesungsfehler des Um- 
kehrpunktes vor, so gehen diese in zwei Bogenmessungen mit gleicher Größe ein. Der Fehler 
kann in diesem Fall bis auf das y 2fache des oben gegebenen Wertes ansteigen. Man muß 
also bestrebt sein, die Ablesefehler klein zu halten. Diese Überlegung gilt auch für alle 
folgenden Fälle. Bei kleiner Dämpfung können wir geradlinigen Amplitudenabfall voraus- 
setzen (Abb. 2). Wir bilden dann das Mittel zwischen drei aufeinanderfolgenden Umkehr- 
punkten, wobei die beiden auf derselben Seite liegenden nur halb bewertet werden. Es 
ergeben sich n =» — 2 Mittelwerte: 


a, a. 
) u 


; USW. 


Daraus erhält man die Gleichgewichtslage: 


I; 
n 
| 
+ a,) + 4, 


v—2 


D. h.: Bei der Mittelbildung ist die erste und letzte Amplitude zu "a, die 
weite und vorletzte zu ”aı zu bewerten; alle übrigen sind voll zu bewerten. 
Die geringe Bewertung des ersten und letzten Umkehrpunktes ist wieder folgerichtig, denn 
man kennt hier weder Schwingungs- noch Dämpfungsstörungen, während sich diese bei den 
mittleren Umkehrpunkten aus dem Schwingungsbild erkennen lassen. Der Wert G, der 
Gl. 3) geht in @,„ der Gl. (1) über, wenn » groß ist. D. h.: Bei einer großen Anzahl von 
beobachteten Umkehrpunkten, die aufeinanderfolgen, kann die Dämpfung für die Berechnung 
der Glei&hgewichtslage vernachlässigt werden. | 


Dagegen kann der mittlere Fehler, den man aus den Gl. im Fall a und b erhält, bei 
Dämpfung erheblich verschieden werden. Als mittlerer Fehler eines Mittels ergibt sieh: 


" ist die Abweichung eines Mittels Z von der berechneten Gleichgewichtslage. Da 7, ent- 
standen ist durch die Mittelung von zwei Werten 1, (siehe Abb. 2), so kann man den Satz 
aufstellen: 


1 
N 
1 
_ 
S 
r 
> 
+J/ + Ze 
m, = — =, 
L l v 3 
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Der mittlere Fehler der Gleiehgewichtslage wird dann: 


2 2:0: 


Man kann dies sofort schließen aus der Gl. (2). 


Abh. 3. Falle; große Dämpfung. 


Im Falle der großen Dämpfung muß man jeden Schwingungsbogen im Dämpfungsver- 
hältnis x teilen und erhält so die einzelnen Gleichgewichtslagen Z, der Bögen (Abb. 3). 
Dabei bestimmt sich 2 aus der Bedingung: 


x sei eine stetige Funktion der Amplitude, es braucht aber keine Konstante zu sein. Doch 
sei es erlaubt, für die Dauer einer Schwingungsperiode mit einem mittleren konstanten 
Wert x zu rechnen. Dann ergibt sich nach der Abb. 3: 


ad, — 4, 

L=a—ı. 
1l+x, 
x,a,+4, 

Ih, 

1+x, 
%, 4, 
‘USW. 


Daraus bereehnet sich die Gleichgewichtslage: 

- 

2 
(5a) 


Für konstantes = vereinfacht sich die Gleichung: 


| | x 


a). 


Es ist hier wieder bezeichnend, daß sich die Dämpfung bei konstantem # nur in der Bewertung 
des ersten und letzten Gliedes bemerkbar macht. Je größer die Dämpfung ist, um so weniger 
wird der erste Umkehrpunkt und um so stärker der letzte Umkehrpunkt bewertet. Die 
Summe der beiden Glieder hat immer das Gewicht 1. Bei allen mittleren Gliedern fällt für 
die Berechnung der Gleiehgewichtslage die Dämpfung vollkommen heraus. @G. geht in den 
Wert @G,„ über, wenn die Dämpfung null wird; denn hierfür ist z=1. Damit wird 


| 
für Dämpfung null. 
I+x 2 
(. unterscheidet sich bei größerem » kaum von @G,„. Man kann also sagen: Formel 1 
zenügt immer zur Berechnung der Gleichgewichtslage, wenn eine größere Anzahl von Um- 
kehrpunkten, die aufeinanderfolgen, beobachtet ist. 


Der mittlere Fehler berechnet sich ganz ähnlich wie im Falle a. Wir haben nur Jetzt 
zwei Anschlußbedingungen und erhalten deshalb als mittleren Fehler von L: 


| 


A; 
N 
/ 7, Ur 
oo, 
An 
a, — a,— a, 
; USW. 
oa, 


4). 


PT- 


3). 


eh 
en 
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Als mittlerer Fehler der Gleichgewichtslage ergibt sich: 


My; 
v? 


Die Fehlergleichung (6) im Falle ce deckt sich mit der Gl. (2) des Falles a, nur daß jetzt eine 
Anschlußbedingung mehr vorhanden ist. Dagegen unterscheidet sie sich von der Fehler- 
gleichung (4) des Falles b um den Wert y2. Dies beruht darauf, daß wir im Falle b das 
Mittel von 2 Bogenmitteln der Fehlerreehnung zugrunde legten und: deshalb für Ivo’ nur den 


halben Wert erhielten. 


Nun sieht man auch, welcher Fehler in der Vermessungskunde gemacht wird, wo die 
Abweichung A jedes Umkehrpunktes gegen eine idealisierte Schwingungskurve der Fehler- 
bereehnung zugrunde gelegt wird, während ich nur die Abweichung » der Mittellagen gegen 
die Gleichgewichtslage berücksichtige. In Abb. 4 und 5 sind die beiden typischen Fehlerfälle 
dargestellt. Bei Abb. 4 wird die Amplitude gestört, während die Gleichgewichtslage erhalten 
bleibt. Bei meinen Gleichungen geht nur die eine gestörte Amplitude mit ihrem Fehler » in 
die Rechnung ein, während bei Zugrundelegung der Umkehrpunkte gerade die fehlerfreien 
Schwingungen mit den größten Abweichungen 4 belastet werden. In Abb. 5 ist die Gleich- 
gewichtslage gestört, während die Amplitude geblieben ist. Man sieht auch hier, daß meine 
Rechnung richtig den Fehler » ergibt, während die Werte 4 streuen. 


Abb. 4. Abb. 5. 


Man muß also die beiden möglichen Fälle, Störung der Amplitude und Störung der 
Gleichgewichtslage, streng auseinanderhalten. Hier interessiert nur die Gleichgewichtslage 
mit ihren Fehlern, und es dürfen nieht Schwankungen der Amplitude und der Dämpfung als 
Fehler der Gleichgewichtslage gewertet werden. 


Algebraisch ist der Beweis für die Richtigkeit meines Rechnungsverfahrens für den 
“all ohne Dämpfung leicht erbracht. Ist 4 der Fehler eines Umkehrpunktes und » der 
Fehler eines Bogenmittels, so gilt: 


So: 
Nach den Regeln der Vermessungskunde bekommt man einen mittleren Fehler der Gleich- 


gewichtslage: 


24° 


Nach Gl. (2) ergibt sich dagegen ein mittlerer Fehler 


Der erste Teil des Zählers unterscheidet sich dadurch, daß die erste und die letzte Be- 
obachtung bei m, nur halbwertig gezählt sind. Bei größerem » macht dies aber nicht viel 


ew. 
ch 

(—3)(r — 1) ) 
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7_----- 
R 
). | 

e 
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aus. Im Zähler steht aber bei m, das Glied 2.14, und dieses ist von ausschlaggebender 
Bedeutung. 


I. Wählen wir in dem Ausdruck I» | alle Glieder positiv, so wird: 
+22. 


Dies ist der Fall, wenn alle I gleiche Vorzeichen haben, also nach derselben Seite gehen. 
Dies heißt, die Gleichgewichtslage schwankt um den Wert 1, während die Amplituden ihre 
gleiche Größe beibehalten, wie dies in Abb. 5 dargestellt ist. Dann wird: 


2, Werden alle Glieder 1» I negativ, so wird: 


m, wird dann sehr klein. Dies ist der Fall, wenn die Gleiehgewichtslage immer dieselbe 
bleibt und nur die Amplituden schwanken, wie dies in Abb. 4 dargestellt ist. Die beiden 
Extremwerte von Fall I und 2 sind natürlich nie zu erreichen, da im Falle 1 zum mindesten 
bei einer Schwingung, nämlieh bei dem Übergang von einer Gleichgewichtslage zur andern, 
die Amplitude schwanken muß und im Fall 2 beim Übergang von einer Amplitudengröße zur 
andern die Gleiechgewichtslage einmal falsch bestimmt wird. 


3. Sind die positiven und negativen Glieder gleichmäßig verteilt, so wird: 
24450. 


Die Gleichgewichtslage und die Amplitude schwanken etwa um dieselben Werte. Es ergibt 
sich dann: 
m, =m,' YV2. 


Aus «diesen Betrachtungen ersieht man, daß die beiden Fehlerreehnungen völlig ver- 
schiedene Werte ergeben. Da nun bei der vorliegenden Aufgabe nur die Gleichgewichtslagen 
und nicht die Größen der Amplituden und ihre Schwankungen interessieren, so ist es rich- 
tiger, nach den neuen Gleichungen zu rechnen und aus den einzelnen Mittellagen den mittleren 
Fehler zu bestimmen. Man darf nieht die Umkehrpunkte der Fehlerrechnung zugrunde legen 
wie dies von der Vermessungskunde vorgeschlagen wird. 


Abh. ©. Abp. 


Zum Schluß sei darauf hingewiesen, daß in manchen Fällen die Gleichgewichtslage nicht 
mit der Mittellage des Schwingungsbogens zusammenfällt®). Ist nämlich die Funktion der 
rückführenden Kraft (7°) unsymmetrisch zu @ (Abb. 6), so wird die eine Amplitude (a,) größer 
als die andere (a,). Die Gleiehgewichtslage (@) ist bestimmt durch die Bedingung P= 0. 
Sie fällt mit der Ruhelage des Systems zusammen. Der Energiesatz verlangt: 


\+P:.de=\-—P-d«. 

D. I. die in der Abb. 6 schraffierten Flächen müssen gleich sein. Man sieht aus dieser Dar- 
stellung, daß die Abweichungen der Bogenmitten (1) von der Gleichgewichtslage mit wach- 
sender Amplitude zunehmen. Trägt man L als Funktion der Amplitude oder der Anzahl n 
der verstricehenen Schwingungen auf, so erhält man in diesem Falle keine um die Gleich- 
gewiehtslage streuenden Messungen, wie es die Wahrscheinliehkeitstheorie verlangt, sondern 
eine Kurve, die sieh mit wachsendem n, d. h. mit abnehmender Amplitude, asymptotisch 
einem bestimmten Werte nähert, der mit der Gleiehgewiechtslage (@) zusammenfällt (Abb. 7). 
Benutzt man die hier angegebenen Reehnungsmethoden für Z, so kann man mit Hilfe einer 
graphischen Darstellung auch in diesen verwickelten Fällen die Gleichgewiehtslage richtig 
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>) Dr, Fromm, Berlin, machte bei der Diskussion des Vortrags hierauf aufınerksam. 
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Das Wirbelsystem der ebenen turbulenten Strömung. 
Von F. Magyar in Wien. 


etrachtet man etwa die von Nikuradse') mittels einer mit der Strömung weitergeführten 
Kamera aufgenommenen Bilder einer turbulenten Flüssigkeitsbewegung, so erscheint zu- 
nächst jeder Versuch, in die Verhältnisse einzudringen, wenig aussichtsreich. Man erkennt 
wohl eine gewisse Symmetrie bezüglich der mittleren Stromschicehte und muß alle Zweifel 
an der Existenz einer der Hauptströmung überlagerten Zusatzströmung endgültig fallen lassen; 
mehr aber läßt sich aus den Bildern unmittelbar kaum entnehmen. Und doch können diese 
eigenartigen, scheinbar regellosen Strömungsvorgänge in gewissem Sinne nur ein Abbild der 
Vorgänge in der Grenzschieht sein, deren Anteil gerade bei der turbulenten Strömung außer- 
ordentlich an Bedeutung gewinnt. | 
Wenn man nun nach Prandtl?’) annehmen kann, daß in der Grenzschicht durch die 
Reibung Wirbel entstehen können, weil in ihr die Geschwindigkeit auf 0 abfällt, so darf man 
nach Lagally*) auch umgekehrt das Wirbelsystem als Ursache der Reibung ansehen. Das 
Wirbelsystem kann dann analog der quellenmäßigen Darstellung als Ursache der ganzen 
Strömung aufgefaßt werden. Die Reibung verringert die Wirbelstärke und bremst dadurch 
die ganze Strömung. Der dem folgenden vorangestellte Grundgedanke, die Untersuchung beim 
Studium des Wirbelsystems zu beginnen, verdankt demnach sein Entstehen in gleicher Weise 
der Prandtlschen Grenzschiehttheorie und der Theorie der Wirbelsehichten von Lagally. 
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Zunächst wurden nun Wirbelverteilungen über @Querschnitte an Hand versuchsmäßig 
bestimmter Geschwindigkeitskurven dureh graphische Differentiation ermittelt. Als Beispiele 
können die Abb. 1 bis 7 dienen. Die Geschwindigkeitskurve der Abb. I ist einer Arbeit von 
Prandtl’) entnommen, während den Abb. 2 bis 7, aufgenommen für verschiedene relative 


Anlauflängen j eines glatten Messingrohres vom Durchmesser d 7,6 em, in welches durch 
einen gut abgerundeten Trichter Luft gesaugt wird, die Arbeit von Schiller und Kirsten’) 


1) J. Nikuradse, Untersuchung über die Geschwindigkeitsverteilung in turbulenten Strömungen. (Göttinger 
Dissertation.) Forsch.A.d.V.D I. Heft 281 und auch Pra ndtl, Nenere Turbulenzforsehung, hydr. Probleme, V.D.1. 
Verlag 1926. 

2) L. Prandtl, Über Flüssigkeitsströmungen bei sehr kleiner Reibung.  Verh. d. TIL. Int. Math. Kongr. in 
Heidelberg 1904. 

3), M. Lagally, Zur Theorie der Wirbelschiehten. Sitz.-Ber. d. kön. bayr. Ak. d. Wiss. math. phyl. Kl. 1015. 

4, Prandtl, Neuere Turbulenzforschung, Hydr. Probleme. V.D.1.-Verlag 1926. 

5) ,. Schiller und H. Kirsten, Die Entwieklung der Gesechwindigkeitsverteilung bei turbulenter Rohr- 
strömung. Z. f. teehn. Phys. 1929 Nr. 7. 
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zugrunde liegt. Mit Ausnahme der für die kleinste Anlauflänge gültigen Abb. 2 erkennt man 
bei allen den linearen Anteil der nahezu laminaren Kernströmung; gegen die Wand zu steigt die 
Wirbelkurve steil an. Allerdings hat man auf diese Weise den Anteil der überlagerten Quer- 
strömung, der das Wirbelbild zu einer Funktion der Längsrichtung und der Zeit macht, noch 
nicht erfaßt. Gerade aber diese der Hauptströmung überlagerte Zusatzbewegung kann allein 
Klarheit über das Wesen «des turbulenten Zustandes bringen. 

Im folgenden wird zunächst eine ebene, nach beiden Seiten ins Unendliche erstreckte 
Strömung zwischen parallelen ebenen Wänden untersucht. Um allen Forderungen gerecht zu 
werden, denkt man sich das Wirbelsystem aus einzelnen Anteilen zusammengesetzt. Die 
innere Flüssigkeitsreibung bedingt die lineare Wirbelkurve der laminaren Strömung, die be- 
kannt und leieht auszuwerten ist. Der Anteil der stark wirbelnden Grenzschicht setzt sich 
zusammen aus einem konstanten Mittelwert der Zirkulation /’ und aus den Schwankungen um 
denselben. Der konstante Mittelwert bedingt eine von Lagally®) beschriebene Parallelströmung 
mit konstanter Geschwindigkeit, die an der Grenzschicht auf den halben Wert absinkt und 
diese mitzunehmen sucht. Die erwähnten Schwankungen der Zirkulation müssen vorhanden 
sein, wenn die immer wieder beobachtete Zusatzbewegung möglich sein soll. Sie treten dort 
auf, wo der Impulsaustausch der Strömung am größten ist, nämlich am inneren Rand der 
Grenzschicht. Die physikalische Realität dieser Annahme wird sofort klar, wenn man über- 
legt, daß eine gleichförmige Grenzschicht als Sonderfall am unwahrscheinlichsten ist. 

'W Setzt man nun mit Rücksicht auf eine Fourier-Entwicklung 
= eine allgemeine periodische Funktion an, so erhält man für die 
TS . Schwankung der Zirkulation: 
TA feos 77; 

Diese Schwankungen haben aber Zusatzgeschwindigkeiten 
v N Y und »’ in der « und y Richtung zur Folge, die ihrerseits wieder 

=?” auf die Grenzschicht zurückwirken’) und die Wirbeldichte an 
| _ f . den Maximalstellen in der idealen Flüssigkeit ständig, in der 
| zähen bis zu einem durch die Reibung bedingten Beharrungs- 
| | wert erhöhen. An den Minimalstellen nimmt umgekehrt die 
Wirbeldichte ab. A, ist ein vorläufig beliebiger Beiwert, e die 


[07 
ES 


Fortpflanzungsgeschwindigkeit der Störungen. Zur Zeit 0 
BE bestimmen sich die Zusatzgeschwindigkeiten mit Hilfe der 
Abb. s. Biot-Savartschen Gleichung wie folgt: 
Substituiert man 
so wird 
In X? + — y)’ | X’+(a’+9y) 
„=> — 
0 
IT 
2n \x: + (a’ + + (a + 
2 \ e 
a+y 


6, Lagallya.a.d. 

s) Die Vorzeiehen der Integrale lassen sieh am besten direkt aus der Abb. S entnehmen. 

»)O0,Schlömileh, Übungsbuch zum Studium der höheren Analysis, Il. Teil, 8. 219, Tenbner 1900, und auch 
Madelung, Die math. Hilfsmittel des Physikers, Springer 1922. 
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infolge des Vorzeichenwechsels von sin n; X im Intervall — & bis 0 verschwinden das zweite 
und vierte Integral der rechten Gleichungsseite und man erhält 


Analog ergibt sich 
+ + = 
\ (Er) cos n, 
E 


x)?’ + (a’ x)? +(a’ + y)? 


Mit der gleichen Substitution wie oben erhält man nach Zerlegung 


+ 
A,T \ Xcosn, X Xsınn,; X 
r 5 
0 
% (a — 4) 
+ & 
X cos n, X  Xsinn,; X 
0 
Somit wird: 
v — A, = (e A, sinn,;,x- Sinn) 


Wie man sich leicht überzeugen kann, genügen diese Zusatzgeschwindigkeiten der Kontinuitäts- 
gleichung 

dur 
Die zusätzliche Bewegung klingt gegen die Mitte zu ab und bleibt am inneren Rande der 
Grenzschicht in der Größenordnung der angenommenen Schwankungen. Wie v. Mises'') 
gezeigt hat, genügen ganz kleine Zusatzgeschwindigkeiten, um das gänzlich veränderte 
Strömungsbild der Turbulenz zu erklären. Die schon wiederholt versuchte Annahme von 
Störungen, die gegen die Wand zu durchweg abklingen, führt auf Widersprüche, da der 
Einfluß der Grenzschicht mit wachsender Entfernung abnehmen muß. In der Grenzschicht 
selbst müssen allerdings die Störungen gegen die Wand zu verschwinden. Diese Frage wird 
im folgenden noch behandelt werden. 


Bei Weglassung der einen Wand ergibt sich aus den Gl. (3) und (4) mit der Substitution 


a+y=z 
direkt 
u 


worin jetzt z den Anstand vom inneren Rand der Grenzschicht der verbliebenen Wand be- 
deutet. Man erkennt, daß die Zusatzgeschwindigkeiten mit wachsendem Abstand von der 
Grenzschicht logarithmisch abklingen. 


Die Abb. 9a und 9b zeigen die Verhältnisse der Kernströmung zur Zeit t=0 in zwei 
um die halbe Periode verschobenen Querschnitten. Um die Hauptströmung, die allen kon- 


stanten Anteilen der Zirkulation entspricht, lagert sich mit wechselndem Zeichen die Zusatz- 


10) Schlömileh.a.a. 
11), R.v. Mises, Elemente der techn. HIydromeehanik, I. Teil S. 71. Teubner 1914. 
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geschwindigkeit «’. Tritt die innere Reibung gegenüber der Grenzschichtwirkung zurück, so 
können auch sattelförmige Geschwindigkeitsbilder gemessen werden, wie sie z. B. Kozeny '"”) 
erwähnt. 

Da die Schwankungen am inneren Rande der Grenzschieht angenommen wurden, läßt 
sich die Strömung in die Grenzschicht analytisch fortsetzen. Die Gl. (3) und (4) lauten hierfür 
entsprechend geändert: 


= ( e a) + e—nk(y+ ar) ) COS 


(6). 


Geht man im Grenzübergang von beiden Seiten an die Schicht y= a’ heran, so findet man 
jene periodisch schwankende Differenz (w — w”), die zur Erklärung der Schwankung_ der 
Zirkulation in dieser Schieht notwendig ist, und außerdem die vom Einfluß der gegenüber- 
liegenden Grenzschieht herrührenden Zusatzgeschwindigkeiten. Die Schicht 7 a’ macht die 
Bewegung mit und die Schwankungen müßten sogar ständig anwachsen, wenn nicht die ge- 
samte schwingende Zusatzbewegung ständig mit Energieverlust verbunden wäre, der nur auf 
Kosten des Druckgefälles hereingebracht werden kann, welches außerdem noch den Reibungs- 
verlust der Hauptströmung zu decken hat. Man kann daher an dem eingangs erwähnten 
Beharrungszustand festhalten, wenn man die Dissipation der Zusatzbewegung vermehrt um 
die Mischverluste gleich jenem Anteil des erfahrungsgemäß zeitlich unveränderlichen Druck- 
gefälles setzt, der nach Abzug des Energiebedarfes der Hauptströmung verbleibt. 


| | | 
Haupr- Haupt- 
stromung 
| 
| | 
| | | 
| 
x Q 
S 
N 
Abb. Abb. vb. Abb. 10, 


Nun haben die Gl. (6) aber noch den Mangel, daß sie «” und ®’” gegen die feste Wand 
zu nieht zum Verschwinden bringen. Diese Schwierigkeit läßt sich aber beseitigen, wenn 
man sich die Auffassung Prandtls'’) zu eigen macht, nach der in der Nähe der Wand nur 
solche Flüssigkeitsballen an der Zusatzbewegung teilnehmen können, deren Durchmesser 
kleiner als der Wandabstand ist. Ein entsprechender Ansatz, der außerdem die Verhältnisse am 


Innenrand der Grenzschieht unverändert läßt, ist z. B., da = proportinal diesem Durch- 


messer Ist, 


)" 

worin d a a’ «die Grenzschichtdieke und m einen zunächst beliebigen Exponenten be- 

deutet, den man in erster Annäherung etwa gleich = oder 1 setzen kann'*). Nun läßt sich 


die Bildung der Schwankungen auch folgendermaßen erklären: sie entstehen aus den gering- 
fügigsten Ursachen an der Wand, breiten sich über die Grenzschicht aus, wobei mit ab- 
nehmendem #7; die Dämpfung kleiner wird '’), und setzen sich am Innenrand der Grenzschicht 


12) ). Kozeny, Über ausgebildete Turbulenz. Sitz.-Ber. d. Ak. d. Wiss. in Wien, Math. naturw. Kl. Abt. Ila. 
37. Bd. Heft 5 u. 6 1928. 

13) ,. Prandtl, Berieht über Untersuchungen zur ausgebildeten Turbulenz. Z. f. ang. M. u. M. Bd. 5. Heft 5. 

14) Verl. Fußnote 18. 

15) Vgl. C.Schaefer, Einführung in die theor. Physik. 8.805. Walter de Gruyter & Co. Berlin-Leipzig 1922. 
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in einem schließlichen Beharrungszustand fest. Praktisch werden jedenfalls überhaupt nur 
iene der ursprünglich übereinander gelagerten Störungen bestehen bleiben, die die größte 
Wellenlänge haben. So wird es erklärlich, daß schon die Annahme einer einfachen cos- 
Funktion in Gl. (1) brauchbare Geschwindigkeitsbilder ergeben kann. 


Für die folgende Überlegung, die sich hauptsächlich auf die Kernströmung bezieht, 
bedeuten: 


!’ die veränderliche Geschwindigkeit der Hauptströmung, 
ihren Mittelwert über den (Querschnitt, 
!/, die zugehörige Geschwindigkeit am Innenrand der Grenzschicht, 


uw — A, Tem 
— (U+w)tl- Sinn,;y 


die Zusatzgeschwindigkeiten und endlich v=U+w, vv die Komponenten der Gesamt- 
geschwindigkeit. 


Da T proportional bez. muß, kann für geschrieben werden, 


wobei jetzt der Beiwert A sowie nach früherem »7,;. mit d von der auf den hydraulischen Radius 
bezogenen Reynoldsschen Zahl R, abhängt. Das Gescehwindigkeitsbild der Hauptströmung 
setzt sich aus Parabeln zusammen, wie etwa die Abb. 10 zeigt. Die zugehörige Wirbelver- 
teilung besteht aus drei Geraden und stimmt mit den aus den Versuchskurven Abb. 2 bis 7 
vermittelten überein, abgesehen von den Abrundungen im Bereich des Innenrandes der Grenz- 
schicht, die aber nach dem vorhergehenden durchaus verständlich sind. Dieses Geschwindig- 
keitsbild dürfte immer gemessen worden sein, weil man annehmen muß, daß die der Zusatz- 
bewegung entsprechenden Vibrationen der Instrumente bei den Versuchen eliminiert wurden. 


Boussinesq hat vorgeschlagen, den Impulsaustausch durch Einführung einer fiktiven 
Schubspannung 

S 
(5) 

zu berücksichtigen, worin der 'Turbulenzkoeffizient e seiner Dimension nach das Produkt 

aus einer Länge mal einer Geschwindigkeit sein muß. Prandtl bestätigt diesen Ansatz'®). 

Setzt man e==1-v’, worin ®’ die (Quergeschwindigkeit und I die Erstreckung der Strömung 

senkrecht zur & — y-Ebene bedeutet, so ist über diese Länge 1 x und y, somit also auch ®’ und vr 

örtlich konstant. Multipliziert man nun die Schubspannung mit 1-d« und integriert über 


die ganze Breite 2a’ der Kernströmung und über die halbe Periode ‚ so ergibt sich für 


den Verlust an Druckenergie über die Länge = 
+ 
—ar 


Dieser Ausdruck ist quadratisch in 7/,. Man findet ihn auch direkt aus den Navier-Stockes- 
schen Gleichungen, ohne ihnen etwas hinzufügen zu müssen. 


Setzt man u—U+w', vv’ und bildet die Ableitungen der rechten Seite, so vereinfacht 
sich die Gleichung; insbesondere wird 


(U+w) 
negativ gleich 


16) Vgl. Fußnote 13. 
11 
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woraus man den gewissermaßen quasistationären Zustand der Strömung erkennen kann. Man 
erhält 


_ 
re 
Da —- selbst veränderlich ist, hat man von dem Mittelwert zu bilden 
+a 
2a n U 
5) 
0 a’ 


Das erste Glied rechts entspricht der reinen Laminarströmung, das zweite den Impulsaustausch 
nach Gl. (9). 


Da die gesamte Zusatzbewegung im Beharrungszustand, also ungedämpft angenommen wurde, 
muß man noch den Anteil der Dissipation 


nK + a’ 
E=x\ | + — 
a’ 
a 
a’ 


der Zusatzbewegung hinzufügen. Damit ist aber auch die Integration über die Länge L 


. . . . JT .. .. 2 
möglich, die sich immer als ganzes Vielfaches von = darstellen läßt. Eine ähnliche Rech- 
k 
nung wäre noch für die Grenzschicht anzustellen. 


Die Bahnen der Flüssigkeitsteilchen infolge der Zusatzbewegung relativ zur Hauptströ- 
mung sind nahezu geschlossen, was man auch aus den Aufnahmen von Nikuradse') trotz 
(der Störung durch die Obertlächenspannung ganz gut entnehmen kann. 


Infolge der starken Dämpfung in ihr bleibt die Strömung in der Grenzschicht fast 
laminar, so daß die Annahme einer parabolischen Geschwindigkeitsverteilung berechtigt er- 
seheint. 10” bleibt dabei, in guter Übereinstimmung mit den gemessenen Werten von 
Niıkuradse'%). Mit abnehmender Reynoldsscher Zahl wird die Grenzschicht immer dieker 
und die Schwankungen an ihrer Oberfläche, die gegen die Mitte zu abnehmen, beeinflussen 
sich immer mehr gegenseitig. Für die Laminarströmung haben sich dann die Verhältnisse 
in der Grenzschicht auf die ganze Strömung ausgebreitet, und die Störungen verschwinden in 
der Mitte durch Interferenz. Demnach stellt sich die laminare Strömung einfach den Grenz- 
fall der allgemeineren turbulenten dar. Es läßt sich ganz elementar zeigen, daß bei der 
Laminarströmung sozusagen die ganze Strömung Grenzschicht ist. Beispielsweise gilt für 
zyl. Rohre 


Ap | 


für die turbulente Strömung 


für die laminare Grenzsehicht.  Bezeichnet U, wieder die Geschwindigkeit am Rande der 
Grenzschicht, so wird an der Stelle y=r 


,. 


17), Vgl. Fußnote 1. 
IS), Kozenya.a. Fig. 6. 
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ın Eliminiert man nun aus den Gl. (13) und (14) das Druckgefälle, so ergibt sich 
3 
1) Ö (D RD (15), 
. 


worin Rp die auf den Durchmesser bezogene Reynoldssche Zahl ist. Vernachlässigt man 
in Gl. (15) ö als klein gegenüber D, so ergibt sich die Gleichung von Kozeny "”) 

80, 43 


2). Diese Vernachlässigung macht aber den Übergang zur Laminarströmung unmöglich, für welche 
nach obiger Überlegung Ö keinesfalls mehr klein ist. Aus Gl. (15) ergibt sich 
32 U, 1 de 
3 
6 Für den Übergang zur Laminarströmung erhält man dureh Gleichsetzen der Druckgefälle 
und mit 
2 4 64 2 
Allgemein ist 
/ x T 
1 Nor, I r 
2 2 2R) U, 
L Ganz analog ergibt sich für die Strömung zwischen parallelen Platten im Abstand 2: 
4a? U, 1 
1- 
wenn N, auf den hydraulischen Radius bezogen ist. Ferner wird 
/ ’r 
öo=a—al/) 
| 0 N, 
t Für den Übergang erhält man wieder 
3x _ . 6 
a 
n und 
e ;. 
n 
A 
u Also wird 
r 
ö=a—-al/ — — = 6. 
Die Grenzschicht wird aber erst in der Anlaufperiode gebildet. Die Wirbelbildung 
beginnt an der Wand und breitet sich gegen die Mitte über jene Grenzschichtdieke aus, die 
der Reynoldsschen Zahl der endgültigen Strömung entspricht. Somit muß die Anlauflänge 
bei der Laminarströmung am größten sein und mit wachsender Reynoldsscher Zahl 
abnehmen. Dies ist tatsächlich beobachtet worden 


Die Schwankungen in der Grenzschieht wurden mit Gl. (1) als allgemeine periodische 
Funktion eingeführt und bieten für eine Annahme noch große Freiheit. Um diese Allge- 
ı meinheit nicht einzuschränken wurde die Rechnung für ein Glied der Fourierschen Reihe 
durchgeführt und die Summierung in den Gl. (9) und (12) über die halbe Periode erstreckt. 
Man wird natürlich von diesem Integrationsintervall unabhängig, wenn man mit Gl. (1) eine 
Reihe einführt, deren Mittelwert über die ganze Periode von O0 verschieden ist. 186 


19) Kozenya.a. ©. 
20), L. Hopf, Zähe Flüssigkeiten. Handb. der Phys. Bd. Vll. S. 140. Berlin, Springer 1927. 
21, Hopf, Handb. der Phys. Band VII, S. 141. 
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Verhalten von Wirbelsystemen. 
Von A. Betz in Göttingen. 
(Aus dem Kaiser Wilhelm-Institut für Strömungsforschung.) 


urch konsequente Anwendung des Kutta-Joukowsky schen Satzes über den Zusammen- 
hang zwischen Tragflügelauftrieb und Zirkulation lassen sich eine Reihe von Sätzen über 
das Verhalten von Wirbelsystemen ableiten. Die Anwendung dieser Gedankengänge auf einige 
Probleme der Tragflügeltheorie gibt Einblick in manche bisher noch wenig beachtete Zu- 
sammenhänge. 
Es sollen nur ebene Strömungsvorgänge behandelt werden. Alle Wirbelfäden sind dem- 
nach gerade und einander parallel (senkreeht zur Strömungsebene). 


l. Allgemeine Sätze. 


1. Der Kutta- Joukowskysche Satz. Wenn ein Körper, um den das Linienintegral der 
Strömung von Null verschieden ist, der also eine Zirkulation 7 besitzt, sich relativ zu dieser 
Flüssigkeit. mit der Geschwindigkeit #» bewegt, so erfährt er eine Kraft senkrecht zur Be- 
wegungsrichtung, die auf die Längeneinheit bezogen 


ist (Kutta-Joukowskvscher Satz, Abb. I). Wenn die Flüssigkeit außer der Zirkulations- 
strömung keine Bewegung ausführt, so ist » die Geschwindigkeit relativ zur Flüssigkeit im 
Unendliehen. Wenn aber die Flüssigkeit außer der Zirkulationsströmung noch andere Be- 
wegungen ausführt, wenn z. B. melirere Wirbel, oder auch Quellen und Senken vorhanden 
sind, so ist nieht ohne weiteres klar, was man als Relativgeschwindigkeit 
F ansehen soll. Es ist nun bekannt, daß für » die Geschwindigkeit des Körpers 
relativ zu jener Strömung zu nehmen ist, welche an der Stelle des Körpers 
bei Abwesenheit desselben herrschen würde. Wir wollen dabei den Körper 
y als unendlich klein voraussetzen, da sonst an verschiedenen Stellen des 
RA1BTZI] Körpers verschiedene Geschwindigkeiten herrschen können. (Dieser Fall kann 
dureh Integration aus unendlich kleinen Körpern aufgebaut werden.) Diese 
schärfere Formulierung des Kutta-Joukowskyschen Satzes kann man 
ohne weiteres einsehen, wenn man weiß, daß ein freier Wirbel, auf den keine 
Kraft wirken kann, sieh mit der gleichen Geschwindigkeit bewegt, welche die Strömung an 
der Stelle des Wirbels bei Abwesenheit des Wirbels haben würde. Der Kutta-Joukowsky- 
sche Satz muß demnach für die Bewegung mit dieser Geschwindigkeit die Kraft Null er- 
geben, d. h. die im Kutta-Joukowsky schen Satz auftretende Geschwindigkeit muß relativ 
zu dieser Bewegung gemessen werden. Man kann dies aber auch durch eine geeignete Ab- 
leitung des Kutta-Joukowskyschen Satzes unmittelbar zeigen. In der Regel wird dieser 
Satz in der Weise abgeleitet, daß man in einem Koordinatensystem, in welchem der Körper 
ruht (stationäre Bewegung), die Drücke berechnet, welche infolge der Überlagerung von Zirku- 
lationsströmung und Translationsströmung auf eine den Körper umschließende Kontrollfläche 
wirken und die Impulse, welche durch sie ein- und austreten. Wenn man nun als Kontroll- 
fläche einen den Körper so eng umschließenden Zylinder wählt, daß die Translationsgeschwin- 
digkeit im ganzen Gebiet der Kontrollfläche als konstant angesehen werden kann, so geht in 
den Kutta-Joukowskvschen Satz als Geschwindigkeit » eben diese Translationsgeschwin- 
digkeit in der Umgebung des Körpers ein, das ist aber die Geschwindigkeit, die an dieser 
Stelle bei Abwesenheit des Körpers herrschen würde. 


Abb. 1. Kutta- 
Joukowskyscher Satz, 


2. Der Schwerpunkt endlicher Wirbelgebiete. Haben wir eine Anzahl von Wirbeln in 
einer Flüssigkeit, so befindet sich jeder einzelne in einer Strömung, welche als Feld aller 
übrigen Wirbel durch deren Größe und Anordnung bestimmt ist, und jeder Wirbel bewegt 
sich mit dieser Strömung. Denken wir uns nun an Stelle der Wirbel einzelne feste Körper 
mit der gleichen Zirkulation wie die Wirbel (etwa rotierende Zylinder), so ist auch die 
Strömung die gleiche. Hindern wir diese Körper an ihrer Bewegung mit der Strömung, 
halten wir sie also an ihrem Orte fest, ohne dabei ihre Zirkulation zu ändern, so erfährt jeder 
Zylinder entsprechend dem Kutta-Joukowskyschen Satz eine Kraft und wir müssen, um 
ihn festzuhalten, eine entgegengesetzte Kraft auf ihn ausüben. Für einen Körper mit der 
Zirkulation 7',, der sich an einer Stelle mit der Geschwindigkeit », befindet, ist diese Kraft 


und steht senkreeht zu ”,. Die Resultierende aus Kräften, die wir auf die Zylinder ausüben, 
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muß an den Grenzen der Flüssigkeit von den Wänden aufgenommen werden, d. h. sie muß 
sleich sein der Resultierenden aus den Drücken der Flüssigkeit auf die begrenzende Wand'). 

Ist nun die Flüssigkeit sehr ausgedehnt, so daß die starr und ruhend angenommenen 
Grenzen derselben überall sehr weit von den Wirbeln bzw. Körpern mit Zirkulation entfernt 
sind, so geht bei fest gehaltenen Wirbeln die Resultierende aus den Drücken auf die Be- 
erenzungswände gegen Null. Die von den Wirbeln herrührenden Störungsgeschwindigkeiten 
nehmen nämlich umgekehrt proportional der ersten Potenz der Entfernung ab, die damit zu- 
sammenhängenden Druckdifferenzen (Bernoullische Gleichung x »?) umgekehrt 


proportional dem Quadrat der Entfernung; die Oberfläche der Begrenzung nimmt linear mit 
der Entfernung zu, so daß sich für die Kraft, die als Produkt aus Druckdifferenz und Fläche 
entsteht, eine Abnahme umgekehrt proportional der Entfernung ergibt. Wenn aber diese 
Kraft auf die Begrenzungswände verschwindet, so muß auch die resultierende Kraft auf unsere 
Körper verschwinden, d. h. es muß sein 
Da die Kräfte P, bzw. die Geschwindigkeiten r, beliebige Richtungen haben können, so ist 
unter 2 eine vektorielle Addition der Kräfte bzw. Geschwindiekeiten zu verstehen. Wir 
können statt dessen auch die Komponenten in der X- und Y-Riehtung gesondert addieren 
und erhalten 
n n n n 
wobei die Zeiger . und y die Komponenten in der X- bzw. Y-Richtung der betreffenden Vek- 
toren bezeichnen. Lassen wir die Körper frei, wobei wir sie wieder dureh gewöhnliche Wirbel 
ersetzen können, so bewegen sie sieh mit den Gesehwindigkeiten v,, und unsere vorstehenden 
Gleichungen bilden eine allgemeine Aussage über die Verschiebung der Wirbel in einer 
Flüssigkeit ohne äußere Kräfte, insbesondere in einer unendlich ausgedehnten Flüssigkeit. Zur 
anschaulichen Formulierung dieser Aussage denken wir uns an Stelle der Wirbel Massen- 
punkte, deren Massen proportional den Wirbelstärken sind. Dabei müssen wir allerdings 
auch negative Massen zulassen, indem die Wirbel mit einem Drehsinn positiven Massen und 
die mit entgegengesetztem Drehsinn negativen Massen entsprechen. Wir können dann von 
einem Schwerpunkt des Wirbelsystemes sprechen, indem wir den Sehwerpunkt des ent- 


fundene Aussage über die Wirbelbewegung folgendermaßen ausdrücken (Satz 1). 

In einer Flüssigkeit, auf welche keine äußeren Kräfte wirken können 
(unendlich ausgedehnte Flüssigkeit), bewegen sich die Wirbel so, daß ihr 
Schwerpunkt relativ zu den starren Flüssigkeitsgrenzen bzw. relativ zu der 
im Unendlichen ruhenden Flüssigkeit unverändert bleibt?). 

Vorausgesetzt ist hierbei natürlich auch, daß außer den betrachteten Wirbeln keine 
weiteren singulären Stellen in der Flüssigkeit vorhanden sind. 

Ist die Flüssigkeit durch starre Wände begrenzt und läßt sich über die resultierende 
Kraft auf die Begrenzungswände. bei festgehaltenen Wirbeln (stationäre Strömung) eine Aus- 
sage machen, so ist damit auf Grund der Gl. (2) eine Aussage über die Bewegung des Wirbel- 
schwerpunktes gegeben. Ist die resultierende Kraft auf die Wände P, so muß 


sein. Bezeichnen wir die relativ zu den starren Wänden mit 
v.,so ist 

1) Durch das Festhalten der Wirbel wird die Strömung stationär (vorausgesetzt, daß außer diesen Wirbeln keine 
weiteren Singularitäten vorhanden und die Begrenzungswände starr und ruhend sind). Man kann daher für die Be- 
reehnung der Drücke die einfache Bernoullische Gleichung p + 4 v2 — konst anwenden. Bei freien Wirbeln ist die 


Form der Strömung in einem bestimmten Zeitpunkt identisch wie bei festgehaltenen Wirbeln, aber sie ist im allge- 
meinen nieht mehr stationär, da die Wirbel ja wandern, also ihre räumliche Anordnung ändern. Damit ändern sich 


aber die Drücke, da bei nieht stationärer Bewegung die verallgemeinerte Bernoullische Gleichung p + 4 v2 +0 1 
— konst gilt ($ = Strömungspotential, bei stationärer Bewegung ist * — 0). Bei freien Wirbeln wird auf die Wirbel 


keine Kraft ausgeübt wie bei festgehaltenen, deshalb muß im seitlichen Durehsehnitt auch die resultierende Kraft auf 
277) 
ot 
und im folgenden benützten äußeren Kräfte auf die Flüssigkeitsgrenzen sind daher nieht die bei freien Wirbeln auf- 
tretenden Drücke maßgebend, sondern jene Drücke, welehe beim Festhalten der Wirbel, also bei stationärer Strömung 
auftreten würden. 

2) Diesen Satz hat bereits Helmholtz in seiner bekannten Arbeit „Uber die Integrale der hydrodynamischen 
Gleiehungen, welehe den Wirbelbewegungen entspreehen“* (Journal f. d. reine und angew. Math. 55, S,25 bis 55, 18) auf 
andere Weise abgeleitet. 


die Begrenzungswände verschwinden. Dies wird dureh die Beschleunigungsdrücke o gerade erreicht. Für die hier 


sprechenden Massensystemes meinen. Unter Anwendung dieses Begriffes können wir die ge- ' 
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Die Schwerpunktsgesch windigkeit ist also 


und senkreeht zur Kraft P gerichtet. Wir erhalten daher den Satz 2: 


Wenn bei festgehaltenen Wirbeln (stationäre Strömung), die von der 
Flüssiekeit auf die Begrenzungswände ausgeübten Drücke eine resultie- 
rende Kraft ergeben, so bewegt sieh der Schwerpunkt der freien Wirbel so, 
wie sichnach dem Kutta-Joukowskyschen Satz ein dort befindlicher Flügel, 
dessen Zirkulation gleich der Summe der Zirkulationen der Wirbel ist, in 
ruhender unendlich ausgedehnter Flüssigkeit bewegen müßte, damit sein 
Auftrieb gleich dieser resultierenden Kraft ist. 


Im allgemeinen werden die Drücke auf die Begrenzungswände und damit ihre resul- 
tierende Kraft nieht ohne weiteres bekannt sein. In manchen Fällen lassen sich aber wenigstens 
gewisse Aussagen leicht machen. Ist z. B. die Flüssigkeit einseitig durch eine ebene Wand 
begrenzt oder zwischen zwei parallele Wände eingeschlossen, so kann die aus den Drücken 
auf die Begrenzungswände resultierende Kraft nur senkrecht zu diesen Begrenzungswänden 
stehen. Da der Schwerpunkt der Wirbel sich senkrecht zu dieser Kraft bewegt, so erhalten 
wir den Satz 3: | 


Befinden sieh Wirbel zwischen zwei ebenen parallelen Wänden, oder 
neben einer einzelnen ebenen Begrenzungswand, so bleibt der Abstand des 
Wirbelschwerpunktes von diesen Wänden unverändert (er bewegt sich pa- 
rallel zu den Wänden)’). 


3. Trägheitsmoment endlicher Wirbelgebiete. Wir denken uns die Wirbel in einer Flüssig- 
keit wieder festgehalten und zerlegen an jedem Wirbel die Geschwindigkeiten in eine radial 
nach dem Schwerpunkt oder von ihm weg gerichtete Komponente und eine dazu senkrechte. 
Ist r der Abstand eines Wirbels mit der Zirkulation I’ vom Schwerpunkt und », die radial 
(nach außen) gerichtete Gesehwindigkeitskomponente, so erfährt dieser Wirbel eine Kraft 


welche senkrecht zu r steht und daher bezogen auf den Schwerpunkt ein Moment T- r bildet. 
Die tangentiale (senkrecht zu r gerichtete) Geschwindigkeitskomponente v; ergibt eine Kraft 
in Richtung von r, welche kein Moment um den Schwerpunkt erzeugt. Die Gesamtheit der 
auf die Wirbel ausgeübten Kräfte kann in eine durch den Schwerpunkt gehende Resultierende 
(aus den von ”, herrührenden radialen Kraftkomponenten) und ein Moment 
M=o2T['v,r 

zerlegt werden. Diese müssen gleieh und entgegengesetzt den auf die Flüssigkeitsgrenzen 
wirkenden Kräften und Momenten sein. Lassen wir die Wirbel frei, so bewegt sich der 


Schwerpunkt entsprechend der Einzelkraft nach den bereits abgeleiteten Gesetzen. Außerdem 
bewegen sich die Wirbel in radialer Riehtung mit den Geschwindigkeiten v,. Da 


or 
ist, erhalten wir 


wobei )/ das Moment der äußeren Kräfte, bei festgehaltenen Wirbeln bezogen auf den Schwer- 
punkt, ist. Ist dieses Null’), so erhalten wir 


I T'r? ist eine Größe, die dem polaren Massenträgheitsmoment Im r? bezogen auf den Schwer- 
punkt entspricht. Wir können sie als Trägheitsmoment des Wirbelsystemes bezeichnen und 
erhalten den Satz 4. 


>) Für eine geringe Anzahl von Wirbeln ist dieses Ergebnis bereits von W. Müller dureh numerische Be- 
rechnung der Wirbelbahnen gefunden worden (Vortrag in der Sitzung der Ortsgruppe Göttingen der Ges. f. angew. 
Math. u. Mech. 10928 und auf der Physikertagung in Prag 1929). 

4, Ob in einem gegebenen Falle ein äußeres Moment herrscht oder nieht, bedarf einer sorgfältigeren Unter- 
suchung als die Frage nach äußeren Kräften, da 2. B. bei Vergrößerung der Grenzflächen zwar die Kräfte mit Iir 
gegen Null abnehmen, aber die Momente endlieh bleiben können, indem bei ihnen der Faktor r als Hebelarm hinzutritt. 


ıl 


ıl 
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Wenn bei festgehaltenen Wirbeln die äußeren auf eine Flüssigkeit 
wirkenden Kräfte kein Moment, bezogen auf denSchwerpunkt des in dieser 
Flüssigkeit befindlichen Wirbelsystems, haben, so bleibt das Trägheits- 
moment dieses Wirbelsystems konstant. 

Wir hätten bei unserer Überlegung an Stelle des Schwerpunktes auch einen beliebigen 
anderen Punkt als Momentenpunkt wählen können. Das Ergebnis wäre aber bei Berück- 
sichtigung des Gesetzes der Schwerpunktsverschiebung das gleiche gewesen. 

Ist das Moment der äußeren Kräfte gleichsinnig mit der als positiv gewählten Dreh- 
richtung der Wirbel, so wächst das Trägheitsmoment gemäß Gl. (6), ist es ungleichsinnig, 
so nimmt das Trägheitsmoment ab. 


4. Wirbelsysteme, deren Gesamtzirkulation Null ist. Die kinetische Energie eines Po- 
tentialwirbels in unendlich ausgedehnter Flüssigkeit in einem Kreisring zwischen r und dr 
und der Schichtdicke 1 ist 5 ( Bradı = . Bei der Integration über die ganze 
\arn a Tr 
Flüssigkeit (von r—=0 bis r = x) ergibt sich sowohl bei Annäherung an r 0 wie an r x 
für die Energie &. Solche Wirbel sind deshalb physikalisch nicht zu verwirklichen. Die 
Schwierigkeit bei r =0 erledigt sich dadurch, daß die physikalischen Wirbel stets einen Kern 


von endlichem Durchmesser haben, in dem die Geschwindigkeit nieht mehr mit 7 gegen x 


zunimmt, sondern endlich bleibt. Die Schwierigkeit bei r = x bleibt aber bestehen (außer 
der Energie ist auch der Drehimpuls &). Man kann deshalb in einer unendlich ausgedehnten 
Flüssigkeit Wirbel immer nur paarweise erzeugen, so daß die Summe der Zikulationen Null 
ist. Das Geschwindigkeitsfeld eines solchen Wirbelpaares nimmt in großer Entfernung um- 
gekehrt wie das Quadrat der Entfernung ab, so daß die Energie auch bei beliebiger Aus- 
dehnung der Flüssigkeit endlich bleibt. Deshalb besteht der Satz 5: 

In einerunendliceh ausgedehnten Flüssigkeit ist dieGesamtzirkulation 
aller Wirbel Null. 


Wirbelsysteme mit endlicher Gesamtzirkulation können nur in endlich begrenzten Flüssig- 
keiten vorkommen. Natürlich kann auch in einer unendlich ausgedehnten Flüssigkeit ein Teil 
der Wirbel in so großer Entfernung sein, daß sie für die Strömung an der betrachteten Stelle 
belanglos sind. Dann können sich dort Wirbelsysteme mit einseitiger Gesamtzirkulation be- 
finden, indem eben die die Gesamtzirkulation zu Null ergänzenden Wirbel sehr weit entfernt 
sind. Da aber Energie und Impuls zweier entgegengesetzter Wirbel mit zunehmender Ent- 
fernung wachsen, so kommen sehr große Entfernungen nur in Fällen sehr großer Energie- 
zufuhr vor. Der Fall, daß man Wirbelsysteme mit der Gesamtzirkulation Null hat, ist des- 
halb verhältnismäßig häufig und bedarf einer besonderen Betrachtung, da der Schwerpunkt 
eines solchen Systemes ja im Unendlichen liegt, so daß sich die vorhergehenden Sätze z. T. 
nicht ohne weiteres anwenden lassen. 

Fassen wir einen Teil der Wirbel zu einer Gruppe zusammen und die übrigen zu einer 
zweiten Gruppe, so können wir jede Gruppe für sich betrachten. Wir können uns der besseren 
Anschauung wegen vorstellen, daß die eine Gruppe die rechtsdrehenden, die andere die links: 
drehenden Wirbel enthält, nötig ist dies aber nicht. Voraussetzung ist nur, daß die Gesamt- 
zirkulation der einen Gruppe gleich und entgegengesetzt der der anderen Gruppe und von 
Null verschieden ist. 

Wenn auf die Flüssigkeit keine äußeren Kräfte und Momente?) wirken, also z. B. bei 
unendlich ausgedehnter Flüssigkeit, müssen die Kräfte und Momente auf die festgehaltenen 
Wirbel Null sein, d.h. die resultierende Kraft auf die eine Gruppe muß gleich und entgegen- 
gesetzt der Kraft auf die andere sein und mit ihr auf einer Linie liegen. Diese Kräfte brauchen 
aber nicht notwendig durch den Schwerpunkt jeder der beiden Gruppen zu gehen. Lassen 
wir die Wirbel frei, so bewegt sich demnach der Schwerpunkt jeder der beiden Gruppen 
senkrecht zu dieser gegenseitigen Kraft und mit gleicher Geschwindigkeit. Wir können dies 
durch den Satz 6 ausdrücken. 

Die Schwerpunkte von zwei Wirbelgruppen mit entgegengesetzt glei- 
cher Gesamtzirkulation bewegen sich parallel zueinander mit gleicher Ge- 
schwindigkeit, behalten also auch konstanten Abstand. 

Da wir zunächst über die Richtung der gegenseitigen Kraft nichts wissen, können wir 
auch über die Richtung der Bewegung nichts aussagen. Wenn die gegenseitige Kraft dureh 
den Schwerpunkt einer Gruppe geht, so ist diese Gruppe frei von äußeren Momenten. Ihr 


5) Bei Wirbelsystemen mit der Gesamtzirkulation Null werden auch die Momente beim Hinausrücken der starren 
Begrenzungswände ohne weiteres beliebig klein (vgl. Fußn. 4 S. 166). 
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Trägheitsmoment ist dann konstant (Satz 4)°). Im allgemeinen wird die gegenseitige Kraft 
nicht genau durch die Schwerpunkte der beiden Gruppen gehen. Wenn aber die beiden 
Gruppen räumlich einigermaßen getrennt und in sich abgeschlossen sind, so liegt die Kraft 
fast immer sehr nahe am Schwerpunkt, und man kann wenigstens näherungsweise die Träg- 
heitsmomente als unveränderlich ansehen. 

Geht die Kraft nieht durch die Schwerpunkte der beiden Gruppen, so ändert sich deren 
Trägheitsmoment. Ist die Kraft aber parallel zur Verbindungslinie der beiden Schwerpunkte 
(S, 5, Abb. 2), was sich darin äußert, daß die Schwerpunkte senkrecht zu ihrer Verbindungs- 
linie wandern, so ist das Moment der Kraft, bezogen auf die beiden Schwerpunkte, entgegen- 
wesetzt gleich. Infolgedessen wächst das Trägheitsmoment der einen Gruppe um ebensoviel, 
wie das Trägheitsmoment der anderen Gruppe abnimmt. (Das eine Trägheitsmoment ist im 
allgemeinen positiv, das andere negativ. Die Absolutwerte der 'Trägheitsmomente nehmen 
dann gleich viel zu oder ab.) Wir können dies in dem folgenden Satz 7 zusammenfassen. 


\ 
/ \ \ 
\ \ 
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Abh. 2. Zwei Wirbeleruppen, deren Schwerpunkte sich Abb. 3. Zwei Wirbelgruppen, deren Schwerpunkte 
senkreeht zu ihrer Verbindungslinie bewegen. sich schräg zu ihrer Verbindungslinie bewegen. 


Bewegen sieh ineiner unendlichausgedehnten, im Unendlichen ruhen- 
den Flüssigkeit die Schwerpunkte zweier Wirbelgruppen entgegengesetzt 
eleicher Gesamtzirkulation senkrecht zur Verbindungslinie der Schwer- 
punkte, so bleibt die algebraische Summe der Trägheitsmomente der 
beiden Gruppen unverändert. 


Bildet die Kraft mit der Verbindungslinie der Schwerpunkte (S, S, Abb. 3) einen Winkel, 
d. h. haben «die Schwerpunkte eine Geschwindigkeitskomponente ”, in Richtung ihrer Ver- 
bindungslinie, so wächst das Trägheitsmoment der einen Gruppe mehr als das der anderen 
abnimmt oder umgekehrt. Auf jeden Fall ändert sieh die Summe der Trägheitsmomente der 
beiden Gruppen. Und zwar ergibt sich für diese Änderung auf Grund der Gl. (6) und der Abb. 3 


Ir?’ bedeutet dabei das Trägheitsmoment der einen Gruppe, das der anderen, 
a ist der Abstand der beiden Schwerpunkte, 27" stellt die Gesamtzirkulation einer Gruppe 
dar. Die Summe der Trägheitsmomente wächst, wenn die Bewegung der Schwerpunkte in 
Richtung von der Gruppe mit positiver Zirkulation nach der Gruppe mit negativer Zirkulation 
hin erfolgt. Bei umgekehrter Bewegung nimmt sie ab. In Gl. (8) ist demnach mit Rücksicht 
auf das Vorzeiehen unter I/’ die Gesamtzirkulation jener Gruppe zu verstehen, welche sich 
nach der anderen hin bewegt. 


II. Anwendungen. 


Bei der praktischen Anwendung dieser Sätze hat man es häufig nicht mit einer Anzahl 
von Einzelwirbeln, sondern mit kontinuierlich verteilten Wirbeln zu tun. Darin liegt aber 
keine Schwierigkeit; es treten einfach [-Ausdrücke an Stelle der 3-Ausdrücke. Etwas unan- 
genehmer ist es, wenn sich die Wirbelsysteme in das Unendliche ausdehnen und dabei viel- 
fach selbst unendlich große Zirkulation aufweisen. Aber auch hierfür lassen sich bei einiger 
Vorsieht die vorstehenden Sätze meist anwenden. 

6), Wir haben zwar bei den früheren Überlegungen und so auch bei Ableitung des Satzes 4 immer vorausgesetzt, 
daß außer dem betrachteten Wirbelsystem keine weiteren Singnmlaritäten, also auch keine weiteren Wirbel vorhanden 
sind. Wenn wir jetzt zwei Gruppen von Wirbeln haben und jede für sieh betrachten, so trifft diese Voraussetzung 
nieht mehr zu. Wir können aber alle früheren Überlegungen in der Weise verallgemeinern, daß wir die Kräfte, die 
zum Festhalten der nieht zur gerade betrachteten Gruppe gehörigen Wirbel nötig sind, als äußere auf die Flüssigkeit 
wirkende Kräfte auffassen. Man überzeugt sich so leieht, daß der Satz 4 in diesem Sinne aueh auf eine Gruppe von 
Wirbeln bei Anwesenheit weiterer Wirbel anwendbar ist. 


\ 
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1. Wirbelsystem hinter einem Tragflüge. Wie die Tragllügeltheorie lehrt‘), bildet sich 
hinter einem Tragflügel bei günstigster Auftriebsverteilung (geringster Widerstand bei ge- 
rebenem Auftrieb) eine Unstetigkeitsfläche, welche an jeder Stelle die gleiche Abwärts- 
veschwindigkeit besitzt. Man kann sich deshalb die Strömung hinter einem Tragflügel so 
vorstellen, wie wenn sich da eine starre Platte, die Unstetigkeitsfläche, mit konstanter Ge- 
schwindigkeit nach unten verschiebt und die Flüssigkeit dabei in Bewegung setzt (Abb. 4). 
Diese Vorstellung ist aber nur in erster Näherung zutreffend, wenn die Störungsgeschwin- 
digkeiten (wozu vor allem die Verschiebegeschwindigkeit w gehört) klein sind gegen die Flug- 
geschwindigkeit. Diese Bewegung wäre nämlich auf die Dauer nur möglich, wenn die Un- 


Abb. 4. Verschiebung einer starren Fläche in einer Abb. >. Aufwieklung der Unstetig 
Flüssigkeit. keitsfläche hinter einem Flügel. 


stetigkeitsfläche tatsächlich starr wäre. Beim Umströmen der Kanten treten seitlich gerichtete 
Saugkräfte P auf die Plattenkanten auf, welche nur von einer starren Platte aufgenommen 
werden können. Bei einer unstarren Unstetigkeitsfläche fallen diese Kräfte fort; die Un- 
stetigkeitsfläche führt infolge dieses Wegfalles der Saugkräfte andere Bewegungen aus, sie 
wiekelt sich von den Rändern ausgehend auf und geht allmählich in zwei getrennte Einzel- 
wirbel über (Abb. 5). 

Ist die Spannweite des Flügels I, so ist nach bekannten Sätzen der Tragflügeltheorie °) 


die Zirkulation je Längeneinheit - für eine solche Unstetigkeitsfläche nach folgender Glei- 
chung über die Spannweite verteilt 

/ . . . . . . . . 
wobei /', die Zirkulation um den Tragflügel in der Symmetrieebene (Mitte) desselben be- 
deutet. Die Abwärtsgeschwindigkeit der Unstetigkeitsfläche vor ihrer Aufwicklung ist 


Wir können die Unstetigkeitsfläche als ein kontinuierlich verteiltes System von Wirbeln mit 
der Gesamtzirkulation Null auffassen. Die Verteilung der Wirbel ist dureh Gl, (9) gegeben. 


Fassen wir die beiden symmetrischen Hälften | =<0 und zu Je einer Gruppe 
zusammen, so ergibt sich aus dem Satz 6, daß der Schwerpunktsabstand der beiden Gruppen 
unverändert bleiben muß. Der Schwerpunkt eines Wirbelsystems nach Gl. (9) von OÖ bis 


liegt, wie man leicht ausrechnen kann, in einem Abstand 


n 1 
von der Mitte. Der Abstand der Schwerpunkte der beiden Gruppen ist demnach 
IT 


Dies ist demnach auch der Abstand der Schwerpunkte der beiden sich bildenden Einzel- 
wirbel (Abb. 5). Im vorliegenden Falle ergibt sich aus der ständigen Symmetrie des Vor- 
ganges, daß die Schwerpunktsverschiebungen stets parallel zur Symmetrieebene stattfinden, 
und daß demnach auch die Einzelwirbel symmetrisch zur ursprünglichen Symmetrieebene 
liegen. 


”), Vgl. etwa Handb. d. Phys. Bd. VII, S. 239 ff. 
s, Handb. d. Phys. a.a.0. 
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Man kann mit ähnlichen Überlegungen auch den zeitlichen Verlauf des Aufwicklungs- | 


vorganges verfolgen, doch erfordert dies ziemlich umständliche Rechnungen. Bis jetzt ist der 
Verlauf des Vorganges in seinen ersten Stadien genauer untersucht, wenn der aufgewickelte 
Teil noch sehr klein gegenüber der ganzen Unstetigkeitsfläche ist?). 


Hier soll versucht werden, über die Größe der Endwirbel und über die Verteilung der | 
Zirkulation in denselben eine näherungsweise Aussage zu gewinnen. Wir teilen die Wirbel £ 
der Unstetigkeitstläche an einer beliebigen Stelle x und fassen die links von dieser Stelle ( 
liegenden Wirbel zu einer Gruppe, die rechts davon liegenden (Abb. 6 diek ausgezogen) zu } 
einer anderen Gruppe zusammen. Wenn wir nun annehmen, daß die gegenseitigen Kräfte 
auf die beiden Gruppen bei festgehaltenen Wirbeln durch den Schwerpunkt der beiden 
Gruppen gehen, was bei der verhältnismäßig starken Konzentration der Wirbel nach den | 


Enden zu und der dadurch bedingten deutlichen Trennung der beiden Wirbelgruppen ziemlich 
genau zutrifft, so muß das Trägheitsmoment einer Wirbelgruppe beim Aufwickeln angenähert 


unverändert bleiben. Die Gesamtzirkulation einer Wirbelgruppe der unaufgewickelten Un- \ 
stetigkeitsfläche von bis Ist 
112 


| 7 
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Abb. 6. Teil der Unstetigkeitsfläche Abbh.7. Zusammenhang zwischen auf- 
und Kreis, auf den sie sich nach dem zewiekelter und nieht aufgewickelter 
Aufwiekeln verteilt. Unstetigkeitsfläche. 


Wir führen für die weitere Rechnung an Stelle der Veränderlichen .« den Winkel 9 ein, der 
nit dureh die Beziehung 


pP = bzw. sing | (7) 
verbunden ist. Es ist also die Wirbelverteilung (Gl. 9) 
öl 
105) 
2 x\? 


Der Schwerpunkt dieser Gruppe hat die Entfernung 


12 
pi \ zde= 2 \ cos dq 4 sing + > (1 
x 


Das Trägheitsmoment der Gruppe in bezug auf die Mitte der Unstetigkeitsfläche (x — 0) ist 
12 


°», MH. Kaden: Aufwieklung einer unstabilen Unstetigkeitsfläche (Göttinger Dissertation). Ing.-Arch. 2 (1931), S. 140. 


\ 
08 
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Das Trägheitsmoment der Gruppe in bezug auf seinen Schwerpunkt ist 


| 


Dieses Trägheitsmoment muß nach dem Aufspulen wieder vorhanden sein. Wir wollen die auf- 
zespulte Gruppe als kreisförmig annehmen, d. h. von der Unsymmetrie absehen, welche durch 
die gegenseitige Einwirkung der beiden aufgespulten Wirbel bedingt ist. Wir können dann 
die Zirkulation als eine reine Funktion des Radius (/"= f(r)) darstellen. Unsere Wirbelpruppe 


von .„ bis , ist zu einer den Kreis mit dem Radius r ausfüllenden Spirale aufgewickelt. 


Es muß dann die Zirkulation 7’, gleich der Zirkulation der ursprünglichen Wirbelgruppe sein 


und außerdem muß das Trägheitsmoment der in diesem Kreis aufgewickelten Wirbel gleich 
dem ursprünglichen Trägheitsmoment der Wirbelgruppe sein 

F 
| 

„2 ( 

Lassen wir r um dr wachsen, so nimmt unter Einhaltung dieser Bedingungen x um d. ab 
bzw. g um dgy zu. Dabei wächst die Zirkulation um 


om, 


und das Trägheitsmoment um 


Nun ergibt sich dureh Differenzieren von Gl. (17) 


0Jx 


| 


2sing 


COS 


1 
+ > sin2g)(l + c0s2y) 


+ > sin2g) 


+4sin?g 
Dies in Gl. (21) eingesetzt, ergibt unter Berücksichtigung von Gl. (20) 


COS 


sın O 
. . . . . . |. 


2r\? l 1 
— pne? - 
(= cos’q + 4sin?g +3 sin24) 


Da sing Y 1 = N ‚ so stellt Gl. (22) den Zusammenhang zwischen r und . dar, d.h. 


sie gibt an, wie groß der Kreis ist, in den ein Stück der ursprünglichen Unstetigkeitsfläche 
übergegangen ist. Da weiterhin die Zirkulation /'—= 7‘, sin g bekannt ist, so ergibt die Gleichung 
auch die Verteilung der Zirkulation im aufgewickelten Endwirbel. In Abb. 7 sind die zu- 
sammengehörigen Werte von I’ und x über r aufgetragen, außerdem die Verteilung der 
2 
Wirbeldichte 
d(r?a) 20T, 
Bildung der betreffenden Grenzübergänge 


Für sehr kleine Werte von 9 ergibt sich aus Gl. (22) durch 


1, 1 
und daraus 
2 
3 (24). 
x 


D. h. ein kleines Randstück der Unstetigkeitsfläche geht bei der Aufwicklung in einen Kreis 


über, dessen Radius 3 der Länge des ursprünglichen Stückes ist. Dieses Ergebnis ist auch 
in der erwähnten Dissertation von Kaden verwandt worden. 
Füry=, ergibt sich >73 Das heißt, der Radius der Endwirbel ist 15: Da 


der Mittelpunkt der Endwirbel im Abstand 75 von der Symmetrieebene liegt, so würde 


das bedeuten, daß sich die beiden Endwirbel gerade berühren. Bei so starker Nähe der 
beiden Wirbel ist aber unsere Annahme, daß die einzelnen Endwirbel kreissymmetrisch sein 


h 
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sollen, sicher nieht mehr zutreffend: zwischen den beiden Wirbeln ist die Geschwindigkeit 
wesentlich höher als außerhalb, dadurch verschieben sich die einzelnen Stromlinien nach 
außen. Die Wirbel werden sich daher in Wirklichkeit nieht berühren. Eine genauere Ver- 
folgung dieser Verhältnisse steht noch aus. Das gefundene Näherungsergebnis mag aber 
wegen seiner Einfachheit ein leicht zu merkendes Bild von der Größenordnung der Wirbel 
geben. Aus Abb. 7 geht hervor, daß der Zusammenhang zwischen r und . ziemlich linear 


2 /I | 
ist. Im größten Teil des Wirbels ist daher gemäß Gl. (24) r- N (3 e), erst im äußeren 
Rand des Wirbels geht der Faktor r in 5 über. Die Kurve für die Verteilung der Wirbel- 


diehte zeigt, daß trotz der verhältnismäßig großen Ausdehnung der Wirbel der Hauptteil der- 
selben doch sehr stark um den Mittelpunkt konzentriert ist. 


2. Vorgänge hinter einem Schaufelgitter. Ähnlich wie bei einem einzelnen Flügel bilden 
sich auch bei einem Schaufelgitter (Abb. 5) hinter den Schaufeln Unstetigkeitsflächen, welche 
bei günstigster Auftriebsverteilung sich relativ zur ungestörten Strömung ebenfalls wie starre 


m 
4 
RA181Z8 IS} 
[RA1s128) = T — 
A | 
Abb. 85. Strömung dureh ein Schaufelgitter mit Abh. 9, Strömung um die theoretischen nieht auf- 
hypothetischen starren Unstetigkeitsflächen. gewickelten Unstetigkeitsflächen im Sehnitt A- A. 


Flächen bewegen würden, wenn die Saugkraft an ihren Rändern aufgenommen werden könnte. 
In Wirklichkeit wiekeln sie sich aber mit der Zeit auf'"). Wir wollen den fast stets vor- 
liegenden Fall betrachten, daß der Abstand a’ der Flächen klein ist gegenüber ihrer Spann- 
weite. Denken wir uns die Unstetigkeitsflächen tatsächlich starr, wie es in der Theorie fast 
immer geschieht, so hätte die Strömung um die starren Flächen weit hinter den Schaufeln 
in der Nähe des Randes den in Abb. 9 dargestellten Verlauf, wenn wir ein Koordinatensystem 
wählen, in welchem diese Flächen ruhen. Da bei diesem Koordinatensystem die Bewegung 
stationär ist, können wir die Bernoullische Gleichung anwenden. Da die Störungsgeschwin- 
digkeit zwischen den Flächen in großer Entfernung vom Rande relativ zu den Flächen ver- 
schwindend klein wird, während außen in der ungestörten Strömung die Relativgeschwindig- 
keit gleich und entgegengesetzt der Verschiebegeschwindigkeit w ist, so haben wir nach der 
Bernoullischen Gleiehung zwischen den Flächen einen Überdruck 


m” 


gegenüber dem Druck in der ungestörten Strömung neben den Flächen ''). Dieser Überdruck 

ist im Gleichgewicht mit den Saugkräften auf die Plattenkanten. Wenn wir nämlich die 

wagerechten Kräfte betrachten, welche auf einen Flüssigkeitsstreifen von der Höhe des 

l"lächenabstandes a’ wirken, so haben wir einerseits infolge der Druckdifferenz zwischen 


n 
außen und innen eine Kraft a’o , ‚, andererseits die Saugkraft auf eine Plattenkante. Impulse 


werden dureh die Grenzflächen nicht transportiert, deshalb muß die Saugkraft 


m? 
sein. 
10) Handb. d. Phys. Bd. VII, S. 272 ff. 
11) Dieht hinter dem Gitter herrsehen andere Drücke und andere Geschwindigkeiten. Dureh Zusammenziehen 


oder Auseinanderrücken der seitlichen Ränder der Unstetigkeitsflächen (positive oder negative Kontraktion) stellt sich 
aber im weiteren Verlauf der Strömung ein Gleiehgewieht mit dem Druck der seitlichen ungestörten Strömung her, 
wobei sich die Geschwindigkeit entsprechend ändert (vgl. Handb. d. Phys. Bd. VII, S. 259). Hier und im folgenden 
betrachten wir nur die Verhältnisse nach diesem Ausgleich. 

Für manche Zweeke ist zu beachten, daß infolge der Breitenänderung der hypothetischen starren Unstetigkeits- 
lächen die Saugkraft auf die Kanten derselben eine Komponente in der Strömungsriehtung besitzt. 


| 
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Bei einem solehen Flächensystem sind die Wirbel sehr stark auf die Ränder konzentriert, 
ın größerer Entfernung vom Rande, also im ganzen mittleren Teil der Flächen sind ja die 
Relativgeschwindigkeiten nahezu Null und damit natürlich auch die Geschwindigkeitsunter- 
schiede auf beiden Seiten der Flächen, die Wirbel. Daher ist die Wirkung der aufgewickelten 
und der nicht aufgewickelten Unstetigkeitsflächen in Entfernungen, die groß gegen den 
Flächenabstand a’ sind, gleich, da ja die räumliche Umlagerung der Wirbel beim Aufwickeln 
unbedeutend ist gegenüber der großen Entfernung. Es ist aber für die Strömung doch ein 
prinzipieller Unterschied, zwischen dem theoretischen Vorgang mit nicht aufgewicklten starren 
Flächen und dem wirklichen mit aufgewickelten Einzelwirbeln, auf den meines Wissens bis- 
her nieht hingewiesen wurde. Er ist im wesentlichen durch den Wegfall der Saugkraft bedingt. 

Die Wirbelgruppe auf der einen Seite befindet sich im Geschwindigkeitsfelde der Wirbel 
der anderen. Wegen der großen Entfernung ändert sich dieses Feld nicht merklich beim 
Aufwiekeln. Wenn wir daher die Wirbel vor und nach dem Aufwickeln uns festgehalten 
denken, so bleibt die gegenseitige Kraft, welche die beiderseitigen Wirbel aufeinander aus- 
üben, unverändert. Die Schwerpunktsgeschwindigkeit der aufgewickelten und der nicht auf- 
gewickelten Wirbel müßte also dieselbe sein. Wenn wir uns aber die Unstetigkeitsflächen 
st: ırr vorstellen, so sind sie nicht mehr kräftefrei, da ja hier noch die Saugkräfte P hinzu- 
kommen. Für die starr gedachten Flächen ist demnach die Geschwindigkeit um einen Betrag 

P 


größer als bei den freien Wirbeln. /’ ist dabei die Zirkulation um den auf einer Seite von 
der Symmetrieebene liegenden Teil einer Unstetigkeitsfläche bezw. um den sich daraus ent- 
wiekelnden Einzelwirbel (im übrigen auch gleich der Zirkulation um die Schaufel in ihrem 
mittleren Teil). Aus Abb. 9 kann man durch Verfolgen des Linienintegrals leicht ablesen, daß 


ist. Damit wird unter Berücksichtigung von Gl. (26) 


Für die Geschwindigkeit der freien Wirbel ergibt sich demnach 


w 


Tatsächlich ist nun auch der Vorgang des Aufwickelns so, daß der Schwerpunkt der an 
jedem Flächenrande angehäuften Wirbel hinter der Geschwindigkeit des mittleren Hauptteiles 
der Flächen zurückbleibt. Während der letztere sich mit der Geschwindigkeit mw bewegt, 


bewegt sich der Schwerpunkt der Wirbel nur mit der Geschwindigkeit , und behält diese 
_ 


Geschwindigkeit auch im Endzustand nach dem Aufwickeln bei. 


Abb. 10. Geschwindigkeiten relativ zur 


ungestörten Strömung. 


TER 
Man kann diese Schwerpunktsgeschwindigkeit , übrigens auch direkt aus dem Felde 


der gegenüberliegenden Wirbel ableiten. Es ist in großer Entfernung identisch mit dem einer 
lotrechten Reihe von konzentrierten Wirbeln, wie nach der Aufwicklung (Abb. 10). Das Feld 
einer solchen Wirbelreihe ist aber in einiger Entfernung eine konstante Geschwindigkeit 


+ + 
richtet ist. Zwischen den beiden Wirbelreihen addieren sich die Felder der beiden Reihen 
zur Geschwindigkeit w, so daß also 


ist. Außerhalb der beiden Reihen haben die Felder entgegengesetztes Vorzeichen, so daß hier 


‚, welche auf der einen Seite nach oben, auf der anderen Seite nach unten ge- 
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die Geschwindigkeit Null ist. Jede Wirbelreihe selbst bewegt sich unter dem Einfluß des 
Feldes der jeweiligen anderen Reihe, also mit der Geschwindigkeit 


m 


Das Aufwickeln hat aber noch eine audere Folge, welche nicht so bequem zu übersehen 


ist wie die Änderung der Wirbelgeschwindigkeit. Bei den starren Flächen hatten wir im 
‚2 


inneren Teil der abgelenkten Strömung einen Überdruck 9=0 ,, dem durch die Saug- 


kräfte an den Flächenkanten das Gleichgewicht gehalten wurde. Nach dem Aufwickeln ist 
keine Saugkraft mehr vorhanden, es kann infolgedessen jetzt auch kein Überdruck im Innern 
zwischen den Wirbelreihen mehr herrschen. Man kann dies auch auf Grund der Bernoulli- 
schen Gleichung einsehen. Bei dem in Abb. 10 gewählten Koordinatensystem ist die Ge- 


n 
schwindigkeit der Innenströmung w, die der Außenströmung Null, die der Wirbel ,„. Um 
stationäre Verhältnisse zu bekommen, müssen wir ein Kordinatensystem wählen, in dem die 


w 
Wirbel ruhen. Dann ist die Geschwindigkeit der Innenströmung _, und die der Außen- 


strömung - (Abb. 11). Beide sind von gleicher absoluter Größe, daher müssen auch die 


Drücke außerhalb und innerhalb der Wirbelreihen gleich sein. 


III 
| 


SIR) .— 


Abh. 11. Gesehwindigkeiten nach der Aufwieklung relativ Abb. 12. Gesechwindigkeitszusammensetzung 
zu den Wirbeln {stationäre Strömung). nach der Aufwieklung. 


Diese Änderung des Druckes während des Aufwickelns ist nun nicht ohne Einfluß auf 
die Strömung im Innengebiet. Wenn wir einen Schnitt durch das Flügelgitter (Abb. S u. 12) 
betrachten und die Bernoullische Gleichung auf die Geschwindigkeiten vor und hinter dem 
Gitter anwenden, so erhalten wir bei Druckgleichheit (aufgewickelte Wirbel, Abb. 12) 


bei einem Überdruck po (nicht aufgewickelte Wirbel) 


Die Geschwindigkeit ist also nach dem Aufwickeln größer als vorher (e,  e,). 

Dieses zunächst wohl etwas merkwürdig erscheinende Ergebnis kaun man sich auch 
noch auf andere Weise klarmachen. Betrachten wir das Gitter mit den abzehenden Wirbeln 
von der Seite, und zwar einmal mit nieht aufgewickelten Wirbelflächen (Abb. S) und einmal 
mit aufgewiekelten Wirbeln (Abb. 12), so ist die Richtung der abgehenden Wirbel verschieden, 


da ihre Eigengeschwindigkeit relativ zur ungestörten Strömung verschieden ist I" und 2]- 
Die Störungsgeschwindigkeit »r, die man Ja als Feld der Wirbel auffassen kann, steht senk- 
recht zu den Wirbeln, hat also in beiden Fällen etwas verschiedene Lage. Bei den nicht 
aufgewiekelten Wirbeln liegen die Wirbel in Richtung von e,, mw steht also senkrecht zu e, 
(Abb. 5) und da e, sich aus der ungestörten Geschwindigkeit e, und der Störungsgeschwindig- 
keit mw zusammensetzt, Ist 

— mw”. 


Im aufgewickelten Zustande bewegen sich die Wirbel mit der Eigengeschwindigkeit =; 


liegen 


daher in ihrer Richtung zwischen e, und e,. Da w und damit auch , wieder senkrecht zu 


den Wirbeln steht, so bilden die Geschwindigkeitsvektoren e,, e,, w ein gleichzeitiges Dreieck, 
dessen Mittellinie die Wirbellinie ist (Abb. 12). Daraus ergibt sich aber, daß e,=e, wird. 181 
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Über die zweckmäßigste Art, 
lineare Gleichungen durch Elimination aufzulösen. 
Von @. Worch in Darmstadt. 


nter dieser Überschrift berichtet Herr Prof. R. Mehmke über das Verfahren von 
B. J. Clasen'!). Dieses Verfahren stellt zweifellos einen originellen und interessanten 
Weg zur Auflösung linearer Gleichungen dar; ob aber für die praktische Rechnung eine 
Ersparnis an Rechenarbeit dabei herauskommt (Herr Mehmke berichtet von einer Ersparnis 
von '/, aller Multiplikationen gegenüber der Gaußschen Methode), erscheint nach dem folgen- 


den wohl noch zweifelhaft. Als Beispiel sei die Lösung von 4 Gleichungen mit den 4 Un- 
bekannten X, bis behandelt: 


Xadan + + den + Aadan = 
Au Ay X. Öec + Xa e= Le 


Das Verfahren von Clasen führt in 3 Gängen zum Ziel. Im ersten Gang werden die 
beiden ersten Gleichungen umgeformt zu 


Xadaaı + + Xadaaı = Zar‘) 
+ Xedenı + Xadanı = Zu: 
Darin bedeuten: 
277 = daa Opa dab 
= — Oda dab 
Zu = 2 Öaa— Zadab: 


Oaaı = DIR 
ba 
Od Öd alba — bı'!ba Aa 
! a 


Zaı — 2a Hba "ba 


Diese beiden Gleichungen mit der dritten ursprünglichen Gleichung werden nun ım 
2. Gang auf die folgende Form gebracht: 


Öu a2 + Aa a2 Zas 


Hierin ist zur Abkürzung gesetzt: 
Öeca dee daaı — Öbe — 
Öde daaı — Iabı — Idaı dac 


He 
Od ha — bı Od c2Veb 
r r r _ . 
Zus = Zvı Hch Der Veh bi 


1, R. Mehmke: Über die zweekmäßigste Art, lineare Gleichungen dureh Elimination aufzulösen. Diese Zeit 
sehrift 1930, H.5, S. 508. 

2) In dieser Form pflegt man in der Baustatik die Elastizitätsgleichungen anzuschreiben. 

3) Nach der Bezeichnungsweise von Mehmke müßte man dar, ea» usw. schreiben. Mit Rücksicht darauf, 
daß sieh diese Abkürzungen als Verschiebungswerte an einem einfach statisch unbestimmten Hauptsystem deuten 
lassen, heiße der dritte Zeiger I (und nicht 2). 
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Ha 
acdaı 
cas 
Zas Hea ea daaı 


Im dritten Gang faßt man schließlich diese drei Gleichungen mit der letzten der ur- 
sprünglichen Gleiehungsgruppe zusammen und formt um zu: 


Ada a3 — 
Zbs 
= 
Aa = 
Hierin ist: 
Ödds = dad — — — 
Zas = 2a — Zee dead — Opa — Zasdad- 


ec» 
r r r Od c2 
eca 
Odds 
Ob ba Hd h Oddz Hd h 
Ohba 
r r Od 
Zus = — ZdsYab 
bs 
aa» 
Das Hda Zds "da > 
ad» 
Die Unbekannten ergeben sich somit zu: 
Oddsz 


Für die praktische Rechnung kann man sich des folgenden Rechenschemas (Zahlentafel 1) 


bedienen. Am linken Rande stehen die Quotienten « und », mit denen — nach der Rechen- 
vorschrift in der 2. Spalte die jeweiligen Zeilen zu erweitern sind. Die zu summierenden 


Zeilen stehen unmittelbar, untereinander; wer ganz sichergehen will, deckt den übrigen Teil 
ab, so daß nur die zu addierenden Werte sichtbar bleiben. 


Wie man leicht nachprüfen kann, sind nach dem Verfahren von Clasen Rechenoperationen 
in folgender Anzahl durehzuführen: 


Quotienten u bzw. 

45 Produkte 

4 Summen mit je 2 Summanden 
3 Summen mit Je 3 Summanden 
2 Summen mit je 4 Summanden 


I (uotienten 

Nieht mitgezählt sind dabei die (in Zahlentafel 1) mit einem * 

sich bei diesen um Rechenproben handelt. 


versehenen Ausdrücke, da es 


Das Verfahren von Gauß: Der Rechnungsgang nach Gauß läßt sich in zwei Abschnitte 
einteilen. Der erste Abschnitt wird ebenfalls in 3 Gängen durchgeführt. Im ersten Gang 
schaltet man die Unbekannte X, aus sämtlichen Gleichungen aus und gewinnt die folgende 
Gleichungsgruppe: 

On hi 7 Aa Od bı 
Xedecı + ei Zei 


da 
| 
PER Aa. 


ur- 


e 


le 


sand 12, Heft 3 
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Zahlentafelil 


> als SQ 
Faktor Rechenvorschrift Vorzahlen la. 
glieder Z 
( 7 Ö ha a a ba 
(2) deh AD 
h (4) = ( 1) h a h a h Aa h 
aa 
7 
(8) = (6) (0) ai Zaı 
ne (11) = de he | — Zur 
(13) = (10) + (11) 
(12) de er Oder bes 
(14) == (5) wen Men” bı 
bı 
Ve} 4 
(15) = (13) ven -ÖderVeb Lea Veh 
hi 
(16) = (14) + (15) Od ba 
He a ) Mea Ögaı Me Mea Bea 
Vena” (18) = — (13) ÖdeaVea Zee Vea 
(19) = (17) +- (18) Ödar 
(22) = — (13) deq Oder Les ded 
Oad (24) = — (19) d I dad dad 
(21)+(22 
+(23)+(24, Od da 
Mic — 5 (26) = (13) ÖecaMdc* 
der 2 
(28) (26) (27) Öeca* Les 
Ödd: 
Ö > r 
(31) = (29) (30) Öbbs” 
!: / r 
da = 32) Zar Mda 
aar 
(33) = (25) va a - Zds Yda 
a2 
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Darin ist zur Abkürzung gesetzt: 


= dep — 

an; = — Öca 

— 

dı d Öya Hd 

Ve dı atd d 


Aus diesen drei Gleichungen entfernt man im zweiten Gang die Unbekannte X,; man 
erhält 
at Nada ca Les 


X. da + da Zds 
worin jetzt bedeuten 


Ob eı 


Ödes bi Hei He, d, 


Zes bei Mei 


dı 


dd: ddı Öpbı 


Zd> Zu Hdı 
Schließlich schaltet man aus diesen beiden Gleichungen noch X. aus und erhält 
Zds 
Darın ist jetzt 
da Od do Od ca Hda 
Der zweite Abschnitt des Gaußschen Verfahrens dient der Ermittlung der Unbekannten 
selbst. Wir brauchen von jedem Gang nur die erste Gleichung. 
Ar An On ar X. + Za 
+ + Aadavı = Zu: 
at Na — Des 


XaOdas Zas 
Daraus folgt: 


Odds 
S, 
X. 
4 4 An (On 
aa 


Wie man erkennt, sind zur Ermittlung der Unbekannten nur die Vorzahlen öy»,, den: , 
und Sowie die Belastungsglieder Z,,, Ze, und Za, nötig. Die übrigen, 
vorhin angeschriebenen Summen sind nur Zwischenwerte, deren besondere Berechnung über- 


| | 
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lüssig ist. Eine Ausnahme machen nur diejenigen Werte, deren Kenntnis zur Berechnung 
ler Quotienten «u erforderlich ist. 

Wir werden praktisch also wie folgt rechnen’): 
Zu: — Za un 
Oaa 
= — — Heı N 
r r r bi 
Ze = Ze — — 
d 
da 
dı 
an — Öca Hd Öebı Hdı Hdı 
r r r da 
Zas = ba — Da Md 
Das Schema für die praktische Durchführung dieser Rechnung zeigt Zahlentafel 2. Der 
Ermittlung der Unbekannten selbst dient die Zahlentafel >. 
Zahlentafel2. 
3jelastungs- 
Faktor Rechenvorsechrift Vorzahlen Be 5 
lieder Z 
(1) du a Öya de a Oda ba 
(2) deh Ödh 
(4)= (2) -+ (3) bi de hi An 
dee Ode be 
Ö 
6)=— — Öyalk de a Me — Öda Me ba 
en ad 
(D)=(5)-+ (6) 
— u; (8) = — (4) e bi Mei Odbı Mei “bi Mei 
01 
(9 (5)+(6)+(8) deca Oder bes 
dad | | 
Na — (11) = (1) ua -Öpa de Hd 
aa 
(12) = (10) (11) 
dı 13 4 r 
(14) — (10) + (11) 
—+- (13) de 
(15) = — (9) | — Zee 
cce2 

(16) = +(13)-+(15) da 

n, 

T- 


4) Vel. u. 


a. Domke: Die Theorie des Eisenbetons. 


Handbneh für Eisenbetonbau. 


4. Aufl., 1. Bd., 8. 522. 


| | 
. 


ISO Woreh. Über die Art. 


lineare Gleichungen durch 


Elimination aufzulösen 
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Math. und Mech. 


Zahlentafel 3. 


ba 
Opa Öya 
bes de hi Xe Ödebı deu 
Odds Oybı Oaa 
Die Abzählung gibt die folgende Anzahl notwendiger Rechenoperationen: 
Zahlen- Zahlen- 
insgesamt 
tafel 2 tafel 3 m 
Quotienten 6 — 6 
Produkte . 20 6 26 
Summen mit je 2 Summanden 6 1 7 
Summen mit je 5 Summanden 4 1 5 
Summen mit je 4 Summanden 2 1 > 
Wuotienten X = 4 4 


In den Elastizitätsgleiehungen in der Baustatik ist nach dem Maxwell-Bettischen 


Satz . 


Das Rechenschema 2 kann 


in diesem Falle noch vereinfachen. 


Wie 


sich leicht nachweisen läßt, ist == die besondere Ermittlung der Zwischenwerte 
On und day, wird also unnötig. 


Zahlentafel 


2a zeigt das Rechenschema für diesen Fall. 


Zahlentafel 2a. 


Faktor Rechenvorschrift Vorzahlen Belastungs 
glieder Z 
(1) Öca Oda Za 
(2) de h 
‘ r 
va (3) : (Dun — Öea Za 
(4) (2) (3) Op bi de bi Od bı AR 
(9) be 
1, (6) (1) u, -Ödalle Me 
Oaa 
der N (4) (4) de br Mei Mei Zpı Meı 
bı 
(9) d 
Hd a ‚10 (1) ua altd La 
hı r 
bı 
Ude — (12) (8) Mda Zea Ude 
Vec» 
4 (9)-+-( 10) 
(13 da 


| 

| 

( 

( 

( 

4 
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Kadner. 


Einfache Konstruktionen zur angenäherten Quadratur Is1 


Die Anzahl der notwendigen Rechenoperationen zeigt die folgende Übersicht: 


Zahlen- Zahlen- 

tafel 2a tafel 3 RER 
(Juotienten 6 6 
Produkte . 16 
Summen mit je 2 Summanden 4 | 5 
Summen mit je 3 Summanden 3 | | 
Summen mit je 4 Summanden > | 5 
(uotienten X = 4 


Beliebig viele Gleichungen. Die Überlegungen, die im vorstehenden für 4 Gleiehungen 
durchgeführt werden, lassen sich unschwer auf den allgemeinen Fall von » Gleichungen aus- 


(lehnen. 


darauf verzichtet. 


Da sich jedoch die allgemeinen Ausdrücke nicht so schnell übersehen lassen, sei 
In der folgenden Übersicht sind nur die jeweils erforderliche Anzahl der 


Reehenoperationen für 5 und 6 Gleichungen zusammengestellt; die Zahlen für 4 Gleichungen 
sind des besseren Vergleichs halber nochmals angeschrieben. 


Übersicht. 
4 Unbekannte 5 Unbekannte - 6 Unbekannte 
Gauß Gaub Gauß 
Clasen Clasen Clasen 

Quotienten . 6 6 14 10 20 15 15 
Produkte. . 45 26 22 84 50 40 140 85 65 
Summen mit je 2 Summanden . 7 ) 25 ) 6 41 11 T 
Summen mit je 3 Summanden { 4 4 7 6 
Summen mit je 4 Summanden 2 3 3 3 5 4 4 7 > 
Summen mit je 5 Summanden 2 3 3 3 iM) ) 
Summen mit je 6 Summanden 2 3 3 
Quotienten X = bzw. 4 4 4 5 5 6 6 6 


Aus dieser Übersicht geht klar hervor, daß das Verfahren von Clasen keine Ersparnis 
an Rechenarbeit gegenüber der Gaußschen Methode mit sieh bringt. Im Gegenteil: Das 
Verfahren von Gauß ist, was den Aufwand an Rechenarbeit anbelangt, dem von Clasen 
vorgeschlagenen Lösungsweg überlegen. 160 


Einfache Konstruktionen zur angenäherten (Juadratur. 


Von F. Kadner in Chemnitz. 


on den Methoden zur angenäherten (Juadratur von Segmenten werden diejenigen mit 

gleichmäßiger Unterteilung des Integrationsintervalles (besonders die Simpsonsche 
Regel) stark bevorzugt. Für die zuerst von Gauß untersuchte Annäherung bei niehtgleich- 
mäßiger Intervallteilung war anscheinend bisher nur theoretisches Interesse vorhanden. 
Daß einfache Sonderfälle dieser Methode sehr wohl praktisch verwendbar sind, soll im folgenden 
kurz gezeigt werden. 


Sind im Integrationsintervall (#„ th) in Punkten „+4; h die Ordinaten 
Y; gezeichnet, und soll eine auf Mittelwertbildung der y; beruhende Flächenformel eine An- 
näherung (p-- 1)-ter Ordnung liefern, d. h. die Fläche einer Parabel (p — 1)-ter Ordnung genau 
angeben, so folgen aus dem Ansatz 
\yd 


| 
en 
1e 
te 
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die bekannten Gleichungen 


+..+6 

| 


für py=2r 
| für py=2r+1 


Aus diesen p Gleichungen mit den 2» Unbekannten 4, €, :. läßt sich 


entweder für py=2n die Flächenformel errechnen, bei der die Annäherung von 
der höchstmöglichen Ordnung p -1l=2n —1 ist, 
oder es lassen sich für y< 2» unendlich viele Flächenformeln zur Annäherung 
(» — 1)-ter Ordnung gewinnen, in denen über 2r p willkürliche Parameter Je 
nach Art der Aufgabe frei verfügt werden kann. 
Bei symmetrischer Intervallteilung A,_;= — Cn-i Sind in (1) die Gleichungen 
mit ungeraden Potenzen der 4; identisch erfüllt, und die verbleibenden Gleichungen 
(2a) für n=2r (2b) für n=2r+1 
1— 


(2a) u. (2b) 


Zee) 


ergeben stets den genauen Wert bis zu einem geraden p, d.h. eine Annäherung von 
ungerader Ordnung. 


Bei höchstmöglicher Annäherung » - 2n kann man auf verschiedene Arten (vgl. Ene. II 
C2, pag. z.B. aus 


1014 

ı 0 4 


die /; und dann aus (1) oder einfacher aus (2a,b) die ec; berechnen. Numerische Lösungen 
hierzu finden sich bei Willers, Methoden der praktischen Analysis, pag. 147. 


Die Lösungen für 


„= md 4, ,=trYD, = 


4 


gestatten folgende einfache Konstruktion: 


\ 


\ 
Abb. 1. Abb. 2 
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je 


en 
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Man teile das Intervall in 3 bzw. 5 gleiche Teile, erriehte über dem mittleren Teil mit den 
Außenteilen als Schenkeln das gleichschenklige Dreieck und trage dessen Höhe von ihrem 
Fußpunkt nach beiden Seiten auf der Grundlinie ab; das ergibt die Punkte x, und «,, und 
die Formel 
F h- 18 


liefert eine Annäherung von 3. bzw. 5. Ordnung. 


Die weiteren Lösungen 


15+2y30 
= 4 = n .) 63 


sind praktisch kaum verwertbar. 


Für geringere Annäherung y  2n mit dem Freiheitsgrad 2n - p ist die allgemeine 
Lösung der Gl. (1) oder (2a)(2b) anscheinend noch nicht versucht worden. 


Der einfachste Fall ist eine Annäherung 3. Ordnung mit einem willkürlichen 
Parameter bei 3 symmetrischen Ordinaten. Die Gl. (2b) lauten hier 


1 C a 
Lund geben aufgelöst 
| 
| Cm DWE 
2 1 
Daraus folgt z. B. für 
3 2 
3 s Em S 
l 1 


Bei Mitverwendung der beiden Grenzordinaten 4,=+#l1l ist 2n— p=2; es 


läßt sich also eine Annäherung von der Ordnung p—1l=2n—-3 erzielen. Der einfachste 
dieser Fälle ist die Simpsonsche Regel n=-3. In den nächsten Fällen ist anzusetzen 


l Cm l l a 


das gibt aufgelöst 


A, 1 4,=t | y2l Am 0 
n=4; 
> 49 64 


Die zugehörigen Konstruktionen sind leicht ausführbar: 


/ 
/ / 
/ 


Abb. 4. 


4 
lech. 
Abb. 3. 
| 
3 
7; 
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Man teile das Intervall in 5 bzw. 7 gleiche Teile und schlage um den einer Intervallgrenze 
benachbarten Teilpunkt den Kreis, der durch die Intervallmitte geht; dieser Kreis schneidet 
auf der einen Grenzordinate die Strecke ab, die von der Intervallmitte nach beiden Seiten 
abgetragen die Fußpunkte der Zwischenordinaten liefert. Die Formeln 

12 | 180 


ergeben dann eine Annäherung 5. bzw. 7. Ordnung. 


Fz»2h 


bzw. Fz23h 


Läßt man in den soeben aufgelösten Gleichungssystemen 4, als Parameter variabel 
unter Streichung der jeweils letzten Gleichung, so erhält man unendlich viele An- 
näherungen 


3. Ordnung: >. Ordnung: 
| —2+122? - 104° 
| 61-4?) | 30(1-7°) 154?(1—4°) Cm 157°(1- 4°) 
Daraus folgt z. B. für ee 11 si 64 
978 12 /4 gleiche 
(3 gleiche Intervalle) | 3 9730 m 
Simpsonsche Regel\ S 1? 
| zählt «doppelt / — 30 C, 30 Cm 30 


Verfügt man im voraus über sämtliche 4 der Ordinatenfußpunkte, so wird im allgemeinen 
2n die erreichbare Annäherung ist dann nur von der Ordnung p— 1=n 1, bei 
symmetrischer Verteilung einer ungeraden Anzahl von Ordinaten jedoch wegen (2b) von der 
Ordnung 


Legt man nur m Ordinaten im voraus fest, so läßt sich mit den Gl. (1) eine An- 
näherung von der Ordnung p - 1=2n—m-— 1 erzielen. Im einfachsten Falle erhält man 
zu einer festen Ordinate /, aus den ersten drei Gl. (1) die Annäherungen 2. Ordnung 
mit 2 Ordinaten 
1 1 FR 


1 1 1 


2.B. 


N 
| 
+ 
2 


— + 5 As - 


Besonders das erste Beispiel dürfte in manchen Fällen praktisch ver- 


wendbar und ausreichend sein: Zieht man bei gegebener Grenz- 
ordinate y, in 3 Intervallbreite Abstand die zweite Ordinate ,, SO I 
gibt die Formel | 
F=z2h-— 
4 
eine Annäherung 2. Ordnung. | 
2 J 
Man erhält z. B. bei graphischer Auswertung von 1 
Abb. 5. 
JT JT 
| 28 | 8 


| 
| 
| 
| | 3 
‘ 9 
| 
| 
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Zu mehrfacher Anwendung bei Unterteilung des Intervalles eignen sich von den an- 
veführten Formeln besonders die mit gleichen Intervallteilen 


| 
- 


und die mit gleichwertigen Ordinaten 


l 
3 | 
Annäherung 
a l 


| 
ww 
.. 


n—5 V 


3. Ordnung 


” 


Annäherung 3. Ordnung 


12 


” 


Annäherung 5. Ordnung 


=] 
. 


Für n—2 ist die Konstruktion in Abh. I angegeben; für n= 3 zeigt sie die einfache 
Abb. 6; für höhere » wird der Vorteil großer Annäherung durch die umständliche Konstruktion 


aufgehoben. 
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KLEINE MITTEILUNGEN 


über die angenäherte numerische Lö- 
sung des Problems der Wärmeleitung. 
In dieser Mitteilunz betrachten wir die Aufgabe der 
Ermittelung der Funktion «, die im Innern eines 
einfach zusammenhängenden. ganz im Endlichen 
liegenden Bereiches der Laplaceschen Differential- 
eleichung Au 0 genügt und auf dem Rande des 
Bereiches der Randbedingung 


h, k und f sind gezebene Funk- 


hu--k 


unterworfen ist. 

tionen, # und x sind die Randwerte der gesuchten 
on 


Funktion und ihrer Ableitung nach der äußeren 
Normale. Dabei soll durchwese Ak=0 und 
h+k-=-0 sein, damit Existenz und Eindeutigkeit 
der Lösung gesichert werden können!t). 


Um diese Randwertaufgabe zu lösen, benutzen 
wir eine Methode, die in bekannten Arbeiten der 
Herren Courant, Friedrichs, Gerschgorin, Lewy. 
L,iehmann. Runge und Wolf entwickelt ist. und die 
Differentialgleichung Au =0 durch eine entspre- 
ehende partielle Differenzengleichung ersetzt: nach- 
folgende, irgendeines (z. B. des Liebmannschen) 
Verfahrens benützende Auswertung der Funktion, 
die jener Differenzengleichung genügt, gibt die an- 
venäherten Werte der gesuchten Funktion für 
die Gitterpunkte, welche im Innern des Bereiches 
lieren. 

Das von uns angewandte Iterationsverfahren ist 
demjenigen von Herrn Wolf?) für die Lösung des 
Neumannschen Problems ähnlich. 

1. Das Auslösungsverfahren. Es sei 
irgendein Bereich gegeben, dessen Rand aus Eeken 
und Seiten des Gitters gebildet wird. Die Seite 
des Gitterquadrats wählen wir zur Längeneinheit. 


Vel. Frank u. Mises, Differentialgleichungen der 


Physik, 2. Aufl. IS. 662. 
2) Diese Zeitschr., Bd. 6 (1926) S. 118-150, 


Wir bezeiehnen mit r die Randwerte der Funk- 
tion ®, die der Differenzengleichung 


und mit die benachbarten Funktions- 
werte auf dem ersten Innenring von Gitterpunkten 
(Abb. 1). 


wo die Steigung d der Ableitung 3 (nach der 
en; 
oo» 
inneren Normale gerichtet) oder — 3, propor- 
N 


tional ist. 


On; 


Die Randbedingung (1) kann man in folgender 
Weise schreiben: 


hr—ıd=!; 
oder 
|) 
(4) 


Wir wollen die dieser Randbedingung zenügende 
Funktion ® durch sukzessive Annäherungen er- 
mitteln. 


Es sei die n-te Näherungsfunktion. r, und 
©, seien entsprechende Werte von r und ww. Dann 


| 
S 
32 
| 
| 
f 
[RA1s316) 
Abb. 6. 
he 
| 
Dann ist 
n 
| | 
| 
| 
A 
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bestimmen wir die Randwerte 7,+,ı von P,+ı in 
fol lt T 


\Mit diesen Randwerten r,„+, lösen wir das 
Dirichletsche Problem, d. h. wir finden ®,+1 
nd bestimmen folglich die w„+,. Dann wird 


Yuto 
"+27 
USW 
I | 
? 
| | | | 
Abb. 1. Abb. °. 


In einspringenden Ecken des Randes definieren 
wir, nach Herrn Wolf3), 


1 
(a+ır! n+1);> 
Wo 
_I+xw, 
"a+ı > 
n 


ist (Abb. 2). 


Wir zeigen nun, daß die Folge r, gegen r kon- 
verziert, wobei r und das entsprechende r der 
Bedinzunz (4) zenüzen. Dann folgt elementar., 
dab die Folge ®, zezen die gesuchte Funktion ® 
mit den Randwerten r konvergiert. 

2. Konvergenzbeweis. Der Beweis wird 
vollbracht, indem wir zeigen, daß im Fall der 
homozenen Randbedinzung (f 


(wo 2 und 9 die zu diesem Fall zehörigen r und 
7 bedeuten) die Anwendung unseres Verfahrens 


ergibt®), und daß die erhaltene Lösung mit den 
Randwerten 2 0 die einzige ist. Daraus folgt 
die Existenz Lösung der Gl. (1) und 
die Konvergenz von r, geren das zugehörige r. 
In der Tat, falls wir festgestellt haben. daß 
? 0 für die homogene Randbedingung (5) die 
einzige Lösung der Gl. (1%) ist, so können wir nach 
dem Hauptsatz der Theorie der linearen Glei- 
chungen behaupten, daß es für die inhomogene 
Randbedingunz eine eindeutire Lösung zibt>). 


h+x' 

5) Die Aufgabe der Auffindung einer Funktion d, welche 
der Gleichung (1’) mit der Randbedingung (5) genügt. ist 
der Lösung des linearen Systems 


ı) In diesem Fall ist ou +14 


äquivalent, wobei die Gleichungen vom Typus (I) den inneren 
Punkten des Bereiches und diejenigen vom Typus (ID) den 
Randpunkten entsprechen. Setzen wir f,k— 0, so erhalten 
wir ein homogenes System, und wenn 2,5 seine einzige 
Lösung ist, so hat das zugehörige inhomogene System aueh 
eine und zwar nur eine Lösung. (Vgl. Courant, Fried 
riehs, Lewy. M. A. Bd. 100, S. 38.) 


Diese Lösung sei mit 2, die zugehörigen Rand- 
werte mit r bezeichnet. Wenn wir dann 


setzen und auf o, unser Verfahren anwenden, 
kommen wir zum Fall der homogenen Bedingung: 


N + (w 4 n) 
r = 
h % 


rw, xg 


eilt dabei (6), so folgt 


Wir müssen also zeigen, daß im Fall der homo- 
eenen Randbedinzung 


h+x 
wir für die Gl. (1”) die einzige Lösung 
P 0 
erhalten. 
Wir setzen -  —=a. Danngibt die Anwendung 
h+x 
unseres Verfahrens 
On 1 . . . . . . . (T). 


Wegen den A und x auferlegten Bedingungen ist 


Den Fall, daß auf dem ganzen Rande #« = 1 (Neu- 
mannsches Problem) ist, erörtern wir hier nicht 
und beschränken uns nur auf Fälle, wo in minde- 
stens einem Randpunkt « < 1 ist. 

Es sei zunächst «< 1 auf dem ganzen Rande. 
Wir setzen 

maxa=A4 (V<A<I) 


max | 9. | = Ay 
max 


Eine der Gl. (1”) genügende Funktion kann ihren 
Maximal- sowie Minimalwert nur am Rande er- 
reichen ®). Deshalb ist 


Ay < Mu . . . . . . . . ( 


Anderseits folgt aus (7) und (8) 


Ontıl=4|Pn 
und 

Weren (9) folgt daher 

und 

AnM,- 

Folglich ist 


M„=P®, 


und desweren auch 


lim 


6, Vel. Liebmann, Münchener Sitzgs.-Ber. 1918. 


| 
| 
| 

lim — r) — 

© 
[17 
| | | 

| 
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Die Eindeutigkeit kann leicht festgestellt wer- 
len, wenn wir die Existenz einer von Null ver- 
:chiedener Lösung © voraussetzen und 


9, 
annehmen. Dann wird offenbar auch 


In 


und (10) wird nicht erfüllt. Für @«e<1 ist der 
Beweis damit erledigt. 

Für den Fall «a<1 bleibt die Beweisanordnung 
bestehen, mit dem Unterschied, daß statt 4 nun 
eine andere Zahl 7, 


A< 


venommen werden muß. 

Für den Neumannschen Fall («a=1 auf dem 
sanzen Rande) kann der Beweis nicht auf die an- 
verebene Weise durchgeführt werden und er- 
fordert Nebenbedingungen (I (w-—r)=0 über 
den Rand). Der für diesen Fall von Herrn Wolf?) 
verebene Beweis für die Konvergenz seines Itera- 
tionsverfahrens setzt die Existenz der Lösung des 
Problems voraus, die Herr Wolf nicht beweist®). 


3. Einnumerisches Beispiel. Als Bei- 
spiel behandeln wir folgende Aufgabe: 

Den Temperaturzustand der Ecke 
einer Ziegelsteinwand bei konstan- 
ter Temperatur des Innenraums 
T=40° © und konstanter Außentem- 
peratur £=09" C zu untersuchen. 


d RA25373 
bb. 3. 


Diese Aufgabe führt zu folgendem Problem der 
Wärmeleitung: Die harmonische Funktion u im 
Bereiche S (Abb. 3) zu bestimmen, mit den Rand- 
bedingungen: 


20.0 u 0,6 


== 

s) Dank einer liebenswürdigen Mitteilung von Herrn 
A. N. Tyehonoff habe ich erfahren, daß er einen Be- 
weis für die Konvergenz des Wolfschen Verfahrens bei 
der Annahme (w—r) gefunden hat, der die Existenz 

der Lösung nieht voraussetzt, und sie folglich beweist. 


für die äußere Fläche (DEF) der Wand und 


 — 300.0 


für die innere Fläche (ABC). 


4 


Abb. 4. 


Um diese Funktion tatsächlich zu berechnen, 
betrachten wir statt des Bereiches S einen neuen 
Bereich 5’ (Abb. 4) mit passend abgeänderten 
Randbedingungen. Diese Randbedingungen lauten: 

75r — 9,6. d = 300,0 (d=w-—r) (11) 
für (ABC) und 


für (DEF). 


Wand 
Razs2 25] 
Abh. 


Auf den Seiten (AD) und (CF) geben wir ein- 
fach die Randwerte r an, die man leicht der Be- 
trachtung des Wärmezustandes der unbegrenzten 
Wand entnehmen kann. Es ist bekannt, daß in 
diesem Falle die Temperatur der Wand geradlinig 
abfällt (Abb. 5): diese Temperaturverteilung ist 
als erste Annäherung bei unserer Bereehnung an- 
venommen. 


Zahlentafell. 


1|32.75132.33 [31.67 [30,50 


28,25]28.24]28.24 27.84] 27,55] 


16.50] 9,24| 7,11 


12,80] 10.16] 9,13] 7,86] 6,47] 5.01| 3,55 


7,65| 7,64| 7,63| 7,60| 7,54] 7.45| 7,30| 7.10| 6,83] 6.45] 5,96| 5.35| 4.63! 3.82] 2,97] 2.11 1,25 


2,47| 2,47] 2,47] 2,46] 2,44] 2.42] 2,36| 2,30] 2,22] 2,09] 1,93] 1,73] 1,50] 1,24! 0,96 | 0,68 | 0,40 


wermun 
serune 
.n.unn 
F 
) 
| 
4 
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Die Zahlentafel I gibt die Resultate der Rech- 
nung®), die nach den aus den Randbedinzungen 
und (12) entsprinzenden Formeln 


1 = 17.55 0,561 w,, 


für (ABC) und 
- 0,324 


für (DEF), durchgeführt worden ist: die Rand- 


werte auf den Seiten (AD) und (CF) blieben wäh- 
rend der ganzen Rechnung unverändert. 


Abb. 6. 


Abb. 7. 


Die Abb. 6 stellt die gewonnene Temperaturver- 
teilung und die Abb. 7 die Isothermen im Innern 


der Wand dar. 
Moskau. D. Panow. 253 


») Alle Rechnungen hat Herr A. Chochloff dureh- 


eeführt. 


On Some Many-Valued Solutions of the 
Equations of Elastic Equilibrium in Polar 
Coordinates. 1. If a body occupying a multiply- 
eonneeted region be dissected along a system of 
barriers which should be introduced to make the 
space oceupied by the body simply-eonneeted and 
if one of the faces of the dissected body formed 
by each barrier be displaced relatively to the other 
by a small rigid-body displacement and then the 
opposite faces joined by removal or insertion of 
thin slices of matter of the same kind as that of 
the body, then the new body so formed will be in 
a state of initial stress. The stress and the strain 
produced in the body must be sinzle-valued and 
eontinuous but the displacements would be dis- 
continuous in cerossing a barrier. That such dis- 
placements can be expressed by many-valued func- 
tions was first pointed out by Weingarten'), whose 
theory was further developed by Volterra and 
Cesaro. Oonsequently any many-valued solution 
of the equations of elastie equilibrium of a body 
oeeupying a multiply-conneected rezion, can be 
applied to the problem of dislocation of that body. 
provided that the stress and strain caleulated from 
such a solution are single-valued and eontinuous 
and the discontinuities of the displacements at 
each barrierrepresent a displacement possible in a 
rigid bodv. In the present paper, many-valued 
solutions have been obtained, of the equations of 
equilibrium in polar eoordinates, and applied to 
problems of dislocation of a portion of a spherical 
shell oceupying a multiplv-conneeted space. 


2. The equations of equilibrium of an elastie 
solid under no bodv forces referred to polar co- 
ordinates, in terms of dilatation and rotation are?) 


u) Sin € —2u 
Sin u 
por 
! S —N 
1 04 
(1) 
1 
0 
Ö. 
(A+-2u)rSin —2u 
S —- () 
in ©) 
where) 
| 
-(2) 
and 
IM). = S () 

2) = sin (8). 
r Sin®\od® 3, Sin } ( 


I, For full referenee on the subjeet see Lo ve' Theory 
of Elastieity (drd. ed.) p. 219. 

Love, loe. eit. p. 139. 

3) Love, loe. eit. p. >6. 
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3. If we assume that up = Pf(r.6) and u,, u 
to be independent of ®, we have 


-) Sin 
br | 
where 


Sin 909% 


Substituting for ©, in (4). we find that z satisfies 
the equation 


182 1 [si 
— — Sin —( 
of whieh a solution is known to be 


| 


where P,„ is Legendre's coeffieient of nth degree 
and u = cos 9, 


Hence 


or f ean be written as 


[cr + Sin (9) 
Therefore 
un —® ( Sin © (10a). 


With the above value of vg and 


1 D 
4, and ug = Sin® P„ . (10b), 


4) ,ove, loe. eit. p. 139. 
5) Love, loe. eit. p. 56. 
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We have also the equations) we have 
V?’(ro,)=0 p 
20, =n(n +1) (Cr"—14 
ıy?(ru,) (/+u)r —2uA=V0 D ap 
4 4 or 206%) - (n+1) Cr" PSin® 
The strain coefficients are given nD 
2 or r 09 1 so that the equations (1) are satisfied. 


We see from (5) that the strain eoeftieients and 
eonsequently the stresses are many-valued unless 
n=1and D=0. 

For single-valued stresses therefore we must 
take 

(11) 
and therefore 
e,r=0, e9a9 = — C cosec? 9, \ 


(12). 
ed ('eosee? = ed r er = 


so that 
everywhere. 


The solution therefore applies to a spherieal 
shell, the spherical surfaces being free from trae- 
tion. Since is infinite when 90 and 
it is obvious that we should exelude the portions 
econtaininze 9=0 and == a from the shell. The 
displacement is many-valued and continuous: it 
can be made one-valued if we make the conven- 
tion 22 +a>PD>a, but it then becomes dis- 
continuous, suddenly deereasinge by the amount 
2aCrSin@ as ® changes from 27 +ato.a. Since 
this change in ««p is proportional to r Sin®, there 
is a dislocation, equivalent to the removal of a thin 
slicee bounded by the planes P=-aFf2x0, and 
subsequent joining of the plane faces so formed. 

If we assume the shell to be bounded by conical 
surfaces and the above displace- 
ment requires for its maintenance, normal trac- 
tions on the conical surfaces of amounts 
— cosee?9, and — cosec? ©, respectively. 

4. For a parallel fissure, let us take the displace- 
ment given bv 


u, in in Sin® | 
A cos 9 cos 
—= AP eos 9Sin . (13 
(13) 
4Sin® 
— > S 
up = A(PcosP + Sin $) Sin?® 


Then we have 


so that the equations (1) are satisfied. Caleulating 
the strain-eoeffieients from (13) we have 


Acos®P 
r Sin? 


4 9) cos 


er e99=— 


r Sin? 
= 4.4 cos “Sin ® 
r Sin? © 
so that 
0. (6). 
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We can therefore apply this solution to a sphe- 
rical shell whose surfaces are free from traection. 
The portions containing 9=V and are 
exceluded, as is obvious from (13), by eonical sur- 
faces and 9—= 9%. The displacement is 
many-valued and continuous, but if we restriet ® 
so that 27 the displacement becomes 
sinzle-valued and discontinuous. Taking only the 
multiply-valued terms in the displacement, we 
have when P=0 


and when P=?2n 


Therefore are eontinuous in erossing the 
barrier at ?=0, while up changes suddenly by 
the amount — 27 A as changes from 27 to 
There is therefore a dislocation due to removal 
of a slice of thiekness 22 at To main- 
tain the displacement given by (13), suitable trac- 
tions must be applied to the conical surfaces 
and 

Suddhodan Ghosh. 127. 


University College of Science. 
Caleutta. 


Kurze Inhaltsangabe: 


Über einige mehrdeutige Lösungen des 
elastischen Gleichgewichtsproblems in 
Polarkoordinaten. 

Nach dem Grundgedanken der Weingartenschen 
Methode werden mehrdeutige Lösungen für ein 
Stück einer Kurelschale abgeleitet, das seitlich 
durch einen Konus begrenzt wird und dessen 

Kugeltlächen spannungsfrei sind. 


Ist Schrägüberschreiten der Fahrbahn 
gefährlich? lie Untersuchung dieser von E. 
Everling in Heft 1, Seite 60/61 dieses Bandes 
aufreworfenen Frage kann dadurch erweitert wer- 
den, daß man das Überschreiten der ganzen 
Fahrbahn in Betracht zieht. Man findet nämlich, 
daß es unter festen zeerebenen Umständen eine 
kleinste Gesehwindigkeit gibt. mit der 
der Fußzänger (oder dgl.) die Straße überqueren 
muß, um noch vor dem herankommenden Fahr- 
zeug auf die andere Seite zu gelangen. und dab 
diese kleinste Geschwindigkeit nur bei einem ganz 
bestimmten Schrägwinkel,. nicht aber beim 
rechtwinkligen Überschreiten möglich ist. 

In den Bezeiehnunzen schließe ich mich an die 
Mitteilung von Everling an. Zu den dort ein- 
veführten kommen noeh hinzu (Abb. 1) die in die 
Straßenrichtung fallende Komponente eb der Ent- 
fernung FW (vom Fußgänger bis zum linken 
Vorderrad des Kraftwagens) und die freie 
Breite des Fahrdamms f-b. In der Abbildung 
stellen die ausgefüllten Rechtecke die Anfangs- 
lage, die schraftierten die Endlage dar. 


Man findet 


cos a + p-b 


und erhält daraus 


I+1-+-peosa 


ecosa—-/sina 
)» ist am kleinsten, wenn 


pl 
esina — feosa 
7+1 


d. h. wenn 


a=arctg —- are sin - (2) 


(+Dye-+P 


ist. 
190° 
W 
RA27EZZ 
Abb. 1. Abb. 2. 


Besonders übersichtlich und bemerkenswert ist 
das Ergebnis, wenn man p=( setzt. Dieser Fall 
ist annähernd verwirklicht, wenn es sich um einen 
den Fahrdamm überschreitenden Fußgänger han- 
delt. Es ist dann der zur kleinsten Geschwindig- 
keit gehörige Schrägwinkel 


Das bedeutet: Ein Fußgänger überschreitet die 
Fahrbahn vor einem herannahenden Fahrzeug 
am bequemsten, d. h. mit möglichst geringer Ge- 
schwindigkeit, wenn er im rechten Winkel 
zu der nach dem linken Vorderrad 
führenden Verbindungsgeradengeht. 
(Abb. 2.) Der zugehörige Geschwindigkeitskoef- 
tizient ist 
+1 


Ein Zahlenbeispiel sei mit f 
a) p=1. b) angefügt: 
a) Als günstigster Schrägwinkel ergibt sich 
aus (2) 
a = 8,53° + 2,130 — 10,66°. 


Der erste Anteil würde dem Überschreiten 
senkrecht zu FW entsprechen, der zweite ist 
der Eintluß von ». Für %k erhalten wir aus 
(1) den Wert 0.247. Bei einer Geschwindig- 
keit des Fahrzeugs von ® = 30 km/h müßte 
demnach der Fahrdamm mit kv = 7,4 km/h 
überschritten werden, 
b) Bei Vernachlässigung von p sind die ent- 

sprechenden Werte 

a, = 8,53° 

Die Mindestgeschwindigkeit ist dann bei 
km/h 


»6km/h. 


e=d, 


Die vorstehende Untersuchung ist davon aus- 
vegangen. daß der Fußgänger zerade noch vor 
dem Kraftwagen vorbeikommen soll. Um zur 
Sicherheit zwischen sich und das Fahrzeug noch 
eine gewisse Strecke 9:b zu brinzen, deren Vor- 
handensein in der Mitteilung von Everling vor- 
ausgesetzt wird, wird er mit einer größeren als 
der Mindestgeschwindierkeit den Fahrdamm iüber- 
qaueren, aber zweckmäßig unter dem günstigsten 
Schrägwinkel, wie er sich nach unserer Unter- 
suchung ergibt. H. Heinrich. 276. 
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BUCHBESPRECHUNGEN 


(Die hier angezeigten Bücher sind durch die VDI-Buchhandlung Berlin NW 7, Ingenieurhaus, zu beziehen.) 


PIETRO BURGATTI, Teoria matematica 
della elasticitä. 3.Bd.der Analisi vettoriale 
venerale e applieazioni. Verlag Nicola Zanichelli, 
Bologna 1951. VIII + 370 >. 


Es handelt sich um den dritten Band eines um- 
fassenden Werkes über die Vektoranalysis, dessen 
ersten Band, Transformazioni lineari Burali-Forti 
und Marcolongo, dessen zweiten Band, Fonda- 
menti di Geometria differenziale Burgatti, Boggio 
und Burali-Forti, Bd. VI, Idrodinamica boggio 
und Bd. V, Elettrieitä e magnetismo Marcolongo 
herausgeben. Vektoranalysis, insbesondere Ten- 
sorrechnung ist das Wurchgängig angewendete 
Hilfsmittel. Sie wird im wesentlichen voraus- 
vesetzt, das einleitende Kapitel von etwa 50 Seiten 
repetiert das Wichtigste aus der Potentialtheorie. 
Die nächsten 60 Seiten sprechen von endlichen und 
unendlich kleinen Verschiebungen. Das folgende 
Kapitel, die Statik, ebenfalls etwa 60 Seiten lang, 
enthält die Statik und stellt den Schwerpunkt des 
Werkes dar. Der stark mathematischen Tendenz des 
Ganzen entsprechend finden sich hier neben der 
Aufstellung der Grundgleichungen die allgemeinen 
Theoreme, so die von Volterra umd Somigliana, «lie 
Methoden von Green und Tedone, sowie (das Rezi- 
prozitätstheorem von Betti. Den Schluß des Kapi- 
tels bildet ein Abriß über Thermoelastizität. Das 
nächste Kapitel bringt dann erst die wichtigsten 
Sonderfälle. Ein besonderes Kapitel ist den Pro- 
blemen von Saint-Venant und Clebsch gewidmet. 
Es folgt die Theorie des Gleichgewichtes von 
dünnen Drähten und Platten. Ein ganz kurzes 
Kapitel über anisotrope Körper sagt von diesen 
wenigstens das Wichtigste, den Schluß bildet ein 
Kapitel über die Dynamik der elastischen Körper 
und schließt mit den Formeln von Kirchhoff. Pois- 
son und dem Prinzip von Huyghens. 

Man findet also das Wichtigste in dem Buch, es 
mangelt erfreulicherweise nicht an Literaturanga- 
ben. und wenn auch die italienische Literatur be- 
sonders beachtet wird, so fehlt doch nicht die be- 
deutende ausländische: die allerwichtigsten Be- 
griffe sind viersprachig gegeben. Stichproben 
lassen das Buch korrekt und lesbar erscheinen. 
So verdient es auch in Deutschland Beachtung. 
Manche mehr mathematische Untersuchung, die es 
enthält, findet sich nieht in der bei uns gebräuch- 
lichen Lehrbuchliteratur. 


Berlin. Hamel. 264 
S. TIMOSHENKO, Prof. der Technischen Me- 
chanik an der Universität Michigan. Schwin- 
sungsprobleme der Technik. Ins 
Deutsche übertragen von Dr I. Malkin., 
New York, und Dr. Elise Helly. Wien. VII 
+ 376 S. mit 183 Abb. im Text. Verlag Julius 
Springer, Berlin 1932. Preis geb. 26 M. 


Die deutsche Ausgabe (dieses Buches, das be- 
reits früher in der Zamm. (Bd. 10. S. 202) eine 
eineehende Würdigung gefunden hat, ist sehr zu 
begrüßen, sowohl vom Standpunkt «des Unter- 
richts aus als sicher auch für den praktisch ar- 
beitenden Ingenieur. Sie konnte noch um einige 
Zusätze und Nachträge ergänzt werden, die 
Fragen des Maschinenausgleichs, Schwingungen 
bei konstanter Reibung, u.a. betreffen. Es dürfte 
in der Tat in der deutschen Literatur über den 
Gegenstand kein Lehrbuch geben, das einerseits 
so stark die prinzipiellen Gesichtspunkte der 
Mechanik ausnützt, andererseits doch überall «die 
Probleme unmittelbar aus der technischen Wirk- 
liehkeit nimmt und dabei dureh die knappe Form 
der Darstellung den Stoff so eingehend behandelt. 


Das Buch, das im Original bereits einen großen 
lrfolg hatte, kann auch in der deutschen Be- 


arbeitung wärmstens empfohlen werden, auch 
sei noch die sprachliche Darstellung hervor- 


gehoben, die nicht das Gefühl der Übersetzung 
aufkommen läßt. F. Noether. 28 


Forschungsarbeiten auf dem Gebiete des Inge- 
nieurwesens, 334. Heft: Untersuchungen an der 
Dieselmaschine: Der Spül- und Ladevorgang bei 
Zweitaktmaschinen. Von KURT NEUMANN. 


Das Arbeitsverfahren raschlaufender Zweitakt- 
Vergasermaschinen. Von OTTO KLÜSENER. 
Berlin 1930, VDI-Verlag, 46 S. m. 72 Abb. Preis 
850 M. 


Die Untersuchungen von Neumann sind wie 
so viele vorhergegangene wärmstens zu begrüßen. 
Seit längerer Zeit sind ja gerade die hier behan- 
delten Fragen besonders brennend, die zu über- 
windenden Schwierigkeiten verzögerten jedoch die 
endgültige Lösung. Nun ist aber eine befriedi- 
gende Übereinstimmung der mit den gewählten 
Annahmen erzielten Ergebnisse mit den Versuchs- 
werten erreicht. Die klare und ins Einzelne 
zsehende Darstellung der bei der Spülung auftreten- 
den Vorgänge, die Ausführung genauer Versuche 
und ihre Erörterung, wobei manche überraschende 
Erscheinung festgestellt werden konnte, die Be- 
rücksichtigung auch geringerer, leicht zu über- 
sehender Einwirkungen, wie etwa die Trägheits- 
wirkung der Spülluft und der Auspuffgase, sind 
besonders hervorzuheben. Die technisch so reiz- 
volle Aufgabe, bei nieht genau bekannten physi- 
kalischen Grundlagen dennoch die Wirklichkeit 
angenähert. zu beherrschen, hat Neumann hier 
vorzüglich behandelt, eine Art zu arbeiten, die 
auch dem Unterzeichneten besonders sympathisch 
ist. Dadurch ist gewiß ein wesentlicher Fort- 
schritt in der Erkenntnis der Erscheinungen in der 
Verbrennungsmaschine erzielt worden. 

Auch die Arbeit von Klüsener ist ganz be- 
sonders lehrreich. Auch hier werden außerordent- 
lich gründliche Berechnungen in glücklicher Weise 
mit Versuchsergebnissen verglichen und gezeigt, 
inwieweit die zugrunde zelerten Annahmen be- 
rechtigt sind. Man gewinnt dabei viele sehr wert- 
volle Einblicke, es wird eine Reihe von Fragen ge- 
löst, allerdings werden auch neue aufgeworfen. 
(Gelegentlich wird deutlich gezeigt, wo sich die 


Vorgänge im wirklichen Betrieb von Laborato- 
riumsversuchen unterscheiden und wo auch die 


sorgfältigsten Messungen Ungenauigkeiten zeigen 
müssen. Es ist bewundernswert, mit welcher 
Gründlichkeit und Kenntnis hier die kleinsten 
Einzelheiten berücksichtigt und durchdacht er- 
scheinen. 

Beide Aufsätze werden gewiß dem engeren 
Fachmanne und der einschlägigen Industrie zum 
Nutzen zereichen. Körner. 257 


Prax am 8. Feber 1932. 


Dr. HANS PETER, Privatdozent an der Univer- 
sität Tübingen, Grenzen der Statistik iin 
der Konjunkturforsehung. Ein Beitrag 
zur Kritik der Wirtschaftsprognose. Veröffent- 
lichungen der Frankfurter Gesellschaft für Kon- 
junkturforschung, herausgeg, von Dr. Eugen Alt- 
sehul. *Heft 5, 73 S. Verlag Kurt Schroeder, 
Berlin 1930. 


Bis vor kurzem herrschte in Deutschland nach 


amerikanischem Muster hohe Konjunktur für Kon- 
junkturforschung. Man glaubte, mit Hilfe dieser 
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Methoden die Wirtschaft beherrschen. die Kon- 
junktur stabilisieren oder mindestens ılie Krisen 
mildern. wenn nieht gar vermeiden zu können. 
Die wirtschaftliche Entwieklung der letzten Jahre 
hat diese Auffassung als Illusion erwiesen. Peter 
untersucht die Ursache dieses Versagens und 
reduziert sie auf das Verhältnis der Statistik zur 
Induktion: die heute herrschenden sozialökono- 
mischen Lehren enthalten keine für diesen Zweck 
ausreichende Theorie, und keine ökonomische 
Statistik kann eine solche fehlende Theorie er- 
setzen. Statistische Gesetze wie in der Physik exi- 
stieren zum mindesten in den herrschenden öko- 
nomischen Lehren nicht. Insbesondere vermag 
niemand die mathematische Natur der konjunktur- 
theoretischen Hauptfunktion. des Trend, zu er- 
kennen, Die ganzen rationalen Funktionen nten 
(irmles können noch so schön als Interpolations- 


formeln passen: sie bieten keinerlei Sicherheit zur 
Extrapolation, 

Die beschreibenden Methoden der Konjunktur- 
forschung umd speziell die Statistik können ohne 
Theorie somit nicht einmal sinnvoll angewenidet 
werden. Auch der Ausweg in die Korrelation ver- 
sagt, denn der Nachweis einer hohen Korrelation 
ist kein Beweis für einen realen. kausalen Zu- 
sammenhang. Somit lassen sich Proenosen auf 
statistische Untersuchungen allein nicht gründen. 
Denn hierfür benötigt man nieht nur die Kenntnis 
der Tatsache, daß eine Abhängigkeit besteht, son- 
dern auch ihre Form. Der Verzieht auf Theorie 
macht die moderne Konjunkturforschung ebenso 
unfruchtbar. wie es seinerzeit die sogenannte 
historische Schule der Nationalökonomie war. 


E. J. Gumbel (Heidelberg). 256 


NACHRICHTEN 


Internationaler Mathematikerkongreß Zürich 1932. 


Wie in Heft 4 «des Jahrgangs 1931 schon be- 
richtet wurde, findet in Zürich vom 4. bis 
12, September ein Internationaler Mathematiker- 
kongreb statt. Aus der Liste der allgemeinen Vor- 
träge seien die folgenden genannt: R. Wavre, 
I/aspeet analytique (du probleme des figures 
planetaires, und S. Bernstein, Sur les liaisons 
entre quantites aleatoires. Von den zehn Sek- 
tionen, in «denen «die einzelnen Vorträge gehalten 
werden, betreffen fünf Teilgebiete der ange- 
wandten Mathematik, nämlich Sektion 4: 
Wahrscheinlichkeitsrecehnung und Versicherungs- 
mathematik. Sektion 5: Astronomie, Sektion 6: 
Mechanik und mathematische Physik und Sek- 
tion 7: Mathematisch-technische Wissenschaften. 
Regre Beteiligung seitens der Mitglieder unserer 
Gesellschaft an den Arbeiten des Kongresses ist 
sehr erwünscht. 


Nutzbarmachung wissenschaftlicher Forschung 
für die Industrie. 

Die Göttinger Vertreter der angewandten Mathe- 

matik umd Mechanik haben unter dem Namen 


„Göttinger  Vermittlungsstelle für  ingenieur- 
wissenschaftliche Untersuchungen" (Göttingen, 


Böttingerstraßbe 8) eine Organisation ins Leben 
eerufen,. die auf breiter Basis unter Heranziehung 
mehrerer Forsehungsstätten wissenschaftlich- 
technische Arbeiten im Interesse der Industrie auf 
Antrag der Interessierten ausführen soll. Sie 
will «damit die in der Industrie vorhandenen For- 
schungsstätten ergänzen und die Durchführung 
von wissenschaftlichen Untersuchungen auch für 
mittlere umd kleinere Betriebe ermöglichen. Göt- 
tingen ist für eine solche Vermittlungsstelle be- 
sonders zreeignet, da hier eine grobe Zahl von 
Instituten mit ihren Mitarbeitern zur Verfügung 
stehen, die für solehe Arbeiten besonders in Be- 
tracht kommen: Institut für angewandte 
Mechanik und Kaiser-Wilhelm-Institut für Strö- 
mungsforschung (Leiter: Prof. Dr. Dr.-Ing. e. h. 
l.. Prandtl): Aerodynamische Versuchsanstalt 
(Leiter: Prof. Dr. A. Betz). 

Bei dieser Gelegenheit sei darauf hinzewiesen, 
daß das unter Leitung von Prof. v, Mises stehende 
Institut für angewandte Mathematik, an der 
Universität Berlin  bestimmungsgemäß auch 


Arbeiten mathematischer Richtung für die Zwecke 
(der Industrie oder anderer öffentlicher und pri- 
vater Unternehmungen auf deren Antrag aus- 
rührt. 


Gesellschaft für angewandte Mathematik 
und Mechanik. 


Ortsgruppe Berlin. 

Für den 27. Mai war ein Vortrag von Herrn 
Priv.-Dozt. Dr. St. Beremann (Berlin) über 
„Mehrdeutige Lösungen bestimmter Strömungs- 
aufgaben" angekündigt. Er mußte wegen Er- 
krankung des Vortragenden verschoben werden. 


Prager Mitglieler. 

Am 21. April sprach Hr. Prof. Dr. W. Müller 
(Prag) über „Laminare Ausbreitungsvorgänge in 
zähen Flüssigkeiten“, 

'Anläßlich seiner Ernennung zum Dr. techn. 
e. h. der Deutschen Technischen Hochschule in 
Prag hielt Hr. Prof. Dr. L. Prandt| (Göttingen) 
am 6. Mai in Prag einen Vortrag über „Neuere 
Ergebnisse der Turbulenzforschung“. 

Am 12. Mai sprach Hr. Prof. Bock (Brünn) 


über „Stabilisierung durch kleine  Sechwin- 
gungen‘. 
Persönliches. 


Hr. Dr. Johannes Piceht (Berlin) habilitierte 
sich an der Technisel.en Hochschule Charlotten- 
burg für das Fach der Optik. 


Hr. Prof. ir, ©. B. Biezeno (Delft) wurde an- 
läßlieh der 300-Jahr-Feier der Universität zu 
Amsterdam von der Philosophischen Fakultät 
zum Dr. h. e. ernannt, 


Hr. Prof. Dr.-Ing. W. Gehler (Dresden) ist 
von der Deutschen Technischen Hochschule in 
Brünn zum Dr. Ing. e. h. ernannt worden, 


Hr. Prof. Dr. Suslow (Odessa) feiert am 
31. Mai das fünfzigjährige Jubiläum seiner wissen- 
schaftlichen und pädagogischen Tätigkeit auf dem 
(Gebiete der Mechanik. 

Hr. Dr.-Ing. Kurt Hohenemser habilitierte sich 


an der Universität Göttingen für das Fach der 
angewandten Mathematik und Mechanik. 
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